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ABSTRACT 

The  design  of  an  experiment  is  given  which  employs  an  inverted  coaxial 
diode  geometry  to  investigate  absorption  of  electromagnetic  radiation  below  and  above 
the  plasma  frequency.  Data  obtained  v^en  a  pure  electron  stream  is  introduced 
into  the  system  are  analyzed,  but  cannot  be  completely  explained.  Difficulties 
encountered  in  Introducing  and  ionizing  a  real  gas  are  discussed. 

The  results  of  an  investigation  concerning  the  form  of  the  electron 
distribution  function  in  such  an  experimental  arrangement  are  reported. 

Finally,  the  process  of  free-free  absorption  is  discussed,  and  expressions 
are  given  for  the  absorption  coefficient  in  the  presence  of  ions  and  neutral  atoms 
or  molecules.  These  are  shown  to  be  essentially  equivalent  to  the  results  of  the 
usual  Lorentzian  theory. 
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I 

INTRODUCTION 

This  report  is  a  continuation  of  a  study  (Goodrich,  et  al  1960)  initiated  two 
years  ago  to  investigate  experimentally  and  theoretically  the  absorption  of  electro¬ 
magnetic  waves  by  a  plasma.  The  paramount  necessity  to  account  fully  for  alter¬ 
native  explanations  of  the  communications  blackouts,  both  re-entry  and  thermo¬ 
nuclear,  was  the  primary  motivation  of  the  work. 

The  evident  lack  of  sufficient  esqperimental  results  to  check  and  verify  the 
existing  plasma  theories  over  an  adequate  range  of  parameters  led  to  a  simple  and 
well-defined  plasma  experiment.  An  inverted  co-axial  diode  was  designed  and 
built  as  an  integral  part  of  a  broadband  co-axial  microwave  transmission  line. 
Microwave  attenuation  data  was  to  be  obtained  for  a  pure  electron  stream  and  for 
an  ionized  gas  case  to  check  and  verify  the  existing  plasma  theories. 

An  independent  theoretical  study  was  carried  out  to  examine  alternative 
and  possibly  new  absorption  mechanisms  of  electromagnetic  waves  in  a  plasma. 
This  led  to  a  consideration  of  the  free-free  absorption  at  frequencies  exceeding 
the  plasma  frequency.  The  free-free  absorption  was  the  only  significant 
mechanism  found  to  be  operative,  or  indeed  even  expected.  A  thorough  quantum 
statistical  formulation  of  this  process  was  initiated  and  completed  successfully. 

From  these  statements,  it  is  clear  that  the  work  reported  here  divides 
naturally  into  three  aspects.  The  first  of  these  is  the  experimental  investigation 


1 


—  THE  UNIVERSITY  OF  MICHIGAN  — - - 

4386-2-F 

which  Includes  the  design  of  the  experiment,  the  results  obtained  in  the  laboratory, 
imd  the  discussion  of  these  results  in  terms  of  the  Lorentzian  plasma.  We  have 
not  been  able  so  far  to  account  for  the  absorption  effects  in  a  pure  electron 
stream  either  in  terms  of  a  simple  equivalent  Lorentzian  model  or  by  some 
other  method.  Due  to  some  experimental  difficulties  the  absorption  data  so  far 
obtained  on  the  ionized  gas  is  inadequate  for  presentation  and  analysis.  The 
second  chapter  is  a  theoretical  effort  to  obtain  the  electron  density  distribution 
functions  in  the  tube  for  an  ionized  gas.  This  function  is  of  cardinal  importance 
in  interpreting  the  experimental  electromagnetic  wave  absorption  effects.  The 
third  chapter  is  a  theoretical  Investigation  of  the  free-free  absorption.  Formulas 
for  the  absorption  coefficients  are  presented  and  discussed.  The  conclusion 
reached,  however,  is  that  this  process  is  simply  an  alternative  description  of  the 
usual  process  of  colllsional  loss,  agreeing  with  the  predictions  of  accepted  theory 
in  their  region  of  joint  validity. 
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A  COAXIAL  LOW^DENSITY  PLASMA  EXPERIMENT 
2. 1  Introduction 

The  scattering  of  electromagnetic  radiation  by  ionized  gases  witii  densities 
below  10  ions  per  cm*^  is  of  general  interest  in  the  communication  field.  Detailed 
experimental  data  in  the  vicinity  of  the  idasma  frequency  is  of  particular  interest 
in  revealing  the  importance  of  possible  interaction  mechanisms  and  the  validity  of 
some  commonly  used  formulas  describing  the  electromagnetic  wave  interaction 
with  plasmas. 

Measurements  of  the  interaction  of  microwaves  with  gas  discharges  have 
been  numerous  since  the  early  experiments  of  Tonks  (1931).  Romell  (1951), 
examined  the  differential  scattering  from  a  3.2cm  diameter  plasma  column 
placed  between  two  yagi  antennas.  Herschberger  (1960)  extended  the  experiments 
of  Romell  (1951)  and  Dattner  (1957).  using  a  0.5cm  diameter  discharge  tube 
placed  across  a  3.0in  x  1.5  in.  rectangular  waveguide,  and  between  open  waveguide 
ends.  The  response  to  variable*  plasma  densities  was  observed  at  a  fixed  signal 
frequency  of  2800  Me.  Golant  et  al  (1961)  have  obtained  data  from  the  positive 
column  of  a  hot  cathode  gas  discharge  in  glass  bulbs  axially  inserted  in  a  cylin¬ 
drical  waveguide. 

Interpretation  of  these  and  similar  e:q>eriments  are  complicated  by  the 
boundary  conditions  around  the  plasma  and  by  the  sensitivity  of  the  results  on  the 
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microwave  structure  near  the  plasma.  Furthermore,  they  are  fixed  frequency 
measurements  in  which  the  functional  dependence  on  the  microwave  frequency 
cannot  be  directly  observed. 

The  present  experiment  was  undertaken  with  three  purposes  in  mind;  first, 
to  produce  an  interaction  in  a  simpler  geometry,  i.  e. ,  one  approaching  the  ideal 
of  a  free  space  plane  wave  Incident  on  a  plane  uniform  ionized  gas  barrier;  second, 
to  observe  the  interaction  as  a  continuous  function  of  signal  frequency,  and  third, 
to  initially  test  the  interpretabilit  of  this  arrangement  with  the  simplest  plasma. 

1.  e. .  a  collision-free  electron  cloud  v4iose  density  and  distribution  should  be 
accurately  predicted.  It  is  believed  that  the  first  two  purposes  have  largely  been 
achieved  by  the  apparatus  described  below.  The  use  of  an  electron  doud.  however, 
has  not  necessarily  produced  an  ideal  plasma,  because  electron-conductor  wall 
collisions  have  replaced  gas  collisions.  Extension  of  the  measurements  to  ionized 
gases  is  a  natural  second  step. 

The  construction  of  the  apparatus  was  finished  during  the  first  quarter. 

The  second  quarter  was  spent  in  shakedown  experiments  and  some  modification 
of  the  cathode  heater  assembly.  At  the  beginning  of  the  third  quarter  a  cathode 
of  14cm  length  (nominal)  was  successfully  activated.  Extensive  reflection  and 
transmission  measurements  were  made  on  this  diode  without  any  gas  present. 

Upon  conclusion  of  these  e?q)eriments.  first  nitrogen  and  then  helium  gas  was 
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introduced  into  the  system.  A  rather  strong  cathode  poisoning  resulted  and  no 
useful  results  were  obtained.  The  purity  of  the  gases  was  questioned,  but  was 
found  to  be  adequate.  AVhen  the  gas  was  introduced  into  the  vacuum  envelope  at  a 
pressure  of  1mm  of  Hg,  it  was  necessary,  of  course,  to  turn  off  the  140  liter  ioni  ¬ 
zation  pump.  Sufficient  amounts  of  'dirty' gas  then  evolved  from  the  system  itself  to 
poison  the  oxide  cathode.  This  was  discouraging.  However,  it  was  discovered 
that  the  adsorption  type  fore-pumps  when  chilled  to  liquid  nitrogen  temperature 
did  not  pump  helium  to  any  significant  extent.  Thus  these  pumps  could  be  used 
to  keep  the  helium  atmosphere  clean  for  a  sufficient  length  of  time  to  conduct  the 
experiment.  By;  this  time  however,  the  14cm  cathode  was  poisoned  beyond 
recovery.  A  new  cathode  was  inserted  for  activation  at  the  beginning  of  the  last 
quarter.  While  activating  this  cathode  an  arc  developed  from  cathode  to  anode, 
and  the  anode  was  punctured  releasing  cooling  water,  so  that  it  was  necessary 
to  replace  the  cathode,  the  anode,  and  one  ceramic  plate.  As  a  result 
the  new  20. 3  cm  cathode  was  activated  only  two  weeks  before  the  end  of 
the  fourth  and  last  quarter.  The  pure  electron  results  were  repeated,  but  there 
was  no  time  left  to  carry  on  the  helium  ionization  experiments. 

The  apparatus  is  described  in  some  detail  in  Section  2.2.  The  measure¬ 
ments  on  the  two  cathodes  and  the  analysis  of  the  data  is  given  in  Section  2. 3. 


THE  UNIVERSITY  OF  MICHIGAN 

4386-2-F 


2. 2  Experimental  Design 

The  experimental  apparatus  is  essentially  a  rigid  50  ohm  transmission 
line  with  a  wide  band  RF  generator  at  one  end  and  suitable  detector  and  load  ter¬ 
mination  at  the  other  end.  A  central  15  in.  section  of  this  line  is  a  removable 
coaxial  diode  which  is  operated  in  a  high  vacuum  by  looping  the  line  throuj^  a  36  in. 
deep  by  18  in.  diameter  vacuum  manifold.  The  diode  produces  an  electron  plasma 
region  with  a  length  (d)  to  annular  width  (rQ-rg^)  ratio  of  about  100.  The  complete 
line  is  shown  in  the  scaled  drawing  of  Figure  2.1.  The  physical  arrangement  of 
the  components  around  die  vacuum  console  is  shown  in  Figure  2. 2.  In  the  diode 
section  the  cathode  forms  the  outer  conductor  and  the  anode  the  inner  conductor 
of  the  coaxial  line,  in  contrast  to  the  usual  arrangement.  Thermionic  emission 
is  produced  by  radiant  heating  of  the  outer  conductor  surface  by  a  tungsten  wire 
heater  surrounded  by  a  platinum  reflector.  Anode  heat  from  the  electron  bombar¬ 
dment  is  carried  away  by  water  flow  inside  the  anode  conductor.  The  water  flow 
is  dead-ended  and  thus  requires  only  one  exit  path  from  the  transmission  line. 

As  shown  in  Figure  2.1,  this  exit  path  is  through  an  adjustable  shorted  quarter- 
wave  line  to  eliminate  reflections. 

The  coaxial  line  was  designed  to  minimize  reflections  and  losses  every¬ 
where  in  the  absence  of  the  ion  cloud.  Except  for  short  external  connecting  cables 
the  line  is  a  rigid  structure,  mostly  of  nickel  tubing,  with  the  central  conductor 
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FIGURE  2.1:  SCHEMATIC  OF  THE  EXPERIMENT 
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LEGEND  FOR  FIGURE  2. 1 

A  Concentric  Coolant  THibe  Entrance  and  Exit 

B  Quarter  Wave  Shorting  Stub 

C  Shorting  Stub  Split  Tee  Block 

O  DC  Isolation  Capacitors 

E  Teflon  Vacuum  Tight  Support  Beads 

F  Support  Bead  Cooling  Coil 

G  Tapered  Transitions  to  Input  and  Output  Coaxial  Cables 

H  Vacuum  Tank 

I  Vacuiun  Tank  Cover  Place 

J  DC  Insulated  Vacuum  Flanges 

K  Diode  Centering  and  Tension  Plate 

L  Tension  Plate  Assembly 

M  Tapered  Transitions  to  Coaxial  Diode 

N  Coaxial  Diode,  Heater  landings  and  Heat  Reflector 

O  Return  Path  Coaxial  Line 

P  Quartz  Support  Bead  Housing 

Q  ’'Vacsorb"  Forepumpe 

R  "Vacion"  High  Vacuiun  Pump 


8 


THE  UNIVERSITY  OF  MICHIGAN 

4386-2 -F 


1 


THE  UNIVERSITY  OF  MICHIGAN 

438 6-2 -F 


supported  by  a  minimum  number  of  insulating  beads.  Gradually  tapered  ends  are 
used  at  the  two  external  terminals  and  at  the  diode  ends.  DC  isolation  of  the 
negative  potential  cathode  is  achieved  by  two  large  series  capacitors  of  hig^ 
dielectric  constant. 

The  long,  smooth  coaxial  geometry  offers  several  key  advantages  to 
the  analysis  of  the  experiment: 

(a)  the  waveguide  can  be  made  to  propagate  only  the  simplest  wave, 
the  transverse  electromagnetic  mode  (TEM), 

(b)  "undean"  boundary  conditions  at  the  extremities  of  the  plasma 
region  become  insignificant  with  a  large  length-to-wldth  ratio, 

<c)  the  waveguide  has  a  very  broad  band  response  which  permits 
variation  of  the  transmitted  RF  frequency,  f,  around  any  fixed 
plasma  frequency  fp.  While  the  reverse  operation  is  also  pos¬ 
sible,  enough  uncertainty  exists  in  the  exact  determination  of 
the  plasma  frequency  that  fp  is  better  used  as  a  parameter  than 
a  variable. 

There  are,  also,  practical  difficulties  in  the  desigpi: 

(a)  the  primary  difficulty  is  the  construction  of  the  long,  narrow 
diode  whose  anode  must  be  coaxially  symmetric  yet  unsupported , 
small  in  diameter  yet  capable  of  high  power  dissipation,  and 
whose  cathode  requires  a  long  internal  oxide  coated  surface, 
not  obtainable  from  commercial  tube  manufacturers, 

<b)  the  cathode  is  a  'one-shot*  device  and  the  diode  is  unavailable 
for  adjustment  once  the  cathode  is  activated  inside  the  vacuum 
system. 


and 


and 
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The  Coaxial  Diode.  The  details  of  the  diode  design  are  shown  in 
Figure  2.3.  The  anode  passes  through  a  teflon  vacuum  seal  above  the  vacuum 
cover  plate,  tapers  to  a  small  radius  in  the  active  cathode  region  and  joins  the 
central  conductor  below  a  second  taper.  The  bottom  tip  is  threaded  into  a  reten¬ 
tion  nut  inside  the  inner  conductor.  Tightening  this  nut  apiflies  the  centering  and 
tensioning  action  of  the  anode  stretcher  assembly  shown  in  Figure  2.1.  The 
cathode  is  a  nickel  tube  IS.  75  in.  long,  with  an  inside  diameter  of  0.356  in.  and  a 
wall  thickness  of  0.012  in. ,  which  makes  a  sliding  fit  into  the  ends  of  the  two 
heavy  nickel  taper  sections.  A  12  in.  longitudinal  tungsten  filament  heater  is 
mounted  on  ceramic  insulators  which  are  positioned  inside  the  cylindrical  platinum 
heat  shield.  A  difference  In  thermal  expansion  of  the  cathode  over  the  anode  of  up 
to  0.15  in.  is  allowed  by  the  sliding  joints,  which,  however,  do  not  impair  electri¬ 
cal  continuity. 

Selection  of  the  outer  conductor  as  the  cathode  was  made  to  optimize  the 
radial  uniformity  of  the  electron  density  distribution.  If  p  is  the  electron  density, 
r^  the  cathode  radius,  and  the  diode  current  per  unit  length,  the  cylindrical 
form  of  Child's  law  for  space  charge  limited  operation  can  be  written  as  an 


expansion  in  r 


(2-1) 
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LEGEND  FOR  FIGURE  2.  3 
(All  dimensions  in  inches) 
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.  790  Inside  Dia. 
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where  the  function  ^  both  for  the  regular  and  the  inverted  diode  was  obtained  by 
Langmuir  and  Bldogett  (1923);  m  and  e  are  the  permittivity  of  free  space 
(MKS),  the  electronic  mass,  and  charge.  The  density  function  is  quite  different 
depending  on  whether  r  is  greater  or  less  than  r^.  Figuio  2. 4  is  a  plot  of  p(r) 
versus  r/r^  and  r^/r  for  the  normal  and  inverted  diodes  respectively.  Both 
functions  go  to  infinity  at  r=rQ  but  the  inverted  diode  has  a  minimum  point  and 
reverses  slope.  Hence,  selection  of  the  anode  radius  such  that  tJt  is  about 
3. 0  or  less,  results  in  a  reasonably  uniform  density  over  all  but  the  outer  radii, 
in  contrast  with  the  normal  case.  The  ratio  used  in  the  present  experiment  is 
2.28.  Since  the  plasma  frequency  depends  on  the  half -power  of  p(r)  we  conclude 
that  the  variation  of  fp  around  its  average  is  less  than  8  percent  over  60  percent 
of  the  plasma  area  as  compared  to  50  percent  variation  over  this  same  area  in 
the  normal  diode.  Furthermore,  the  influence  of  the  spike  in  p(r)  near  the 
cathode  is  less  important,  in  the  inverted  case,  in  determining  the  effective 
average  dielectric  constant  for  the  plasma. 

The  radial  dimensions  should  be  as  large  as  possible  for  adequate  coolant 
area  and  for  precision  of  aimular  spacing.  They  are,  however,  severely  limited 
by  power  requirements.  The  companion  equation  to  (2-1)  relating  voltage  and 
charge  density  is 

V(r)=--^  7r2j32  .  (2-2) 

e  4 
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Hence  the  plate  power  per  unit  length  is,  from  (2-1)  and  (2-2), 

The  ratio  rj^/r^  is  approximately  fixed  by  the  need  for  density  uniformity, 
and  the  value  of  p  at  r=rj^  is  typical  of  the  annular  region  average.  Hence,  (2-3) 
displays  the  strong  dependence  of  power  on  r^  and  p  average.  There  is  evidently 
a  need  to  minimize  rg^  to  reach  the  highest  possible  p  .  The  first  diode  was  con¬ 
structed  with  the  values, 

*‘c/*‘a  "  2.28,  and  rg=  0.078  in. 

A  safe  ui^r  bound  for  heat  flux  across  the  anode  wall  is  about  100  watta/cm^.  For 
this  limit  and  the  above  radii,  (2-3)  gives  the  electron  density  at  the  anode 
p(r^  « 1 .6  X  10^/cm^. 

g 

It  is  not  reasonable  to  push  the  electron  density  much  above  2 .0x10  since 
this  would  require  the  radii  to  be  unusably  small  and  high  pressure  cooling  would  be 
necessary.  Electron  densities  near  or  above  10^^/cn:P  can  only  be  reached  by  pulsed 
operation  of  the  diode  plate  voltage.  This  has  not  been  attempted  as  yet. 

The  axial  length  of  the  diode  must  be  long  to  provide  enough  plasma  length 
to  produce  clearly  measurable  changes  in  the  reflection  and  transmission  coef¬ 
ficients,  r  and  T,  in  the  vicinity  of  the  plasma  frequency.  The  anode  cannot  be 
bead  supported  in  and  near  the  active  cathode  region,  however,  because  of 
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sputtering  from  the  cathode.  Hence  the  anode  length  between  closest  bead  sup¬ 
ports  should  be  small  for  mechanical  rigidity.  A  compromise  Is  therefore  neces¬ 
sary  In  the  choice  of  length.  The  effect  of  plasma  length,  d,  on  T  and  can  be 
demonstrated  by  considering  the  simplest  possible  model  of  Interaction,  namely, 
the  passage  of  the  principal  coaxial  mode  (TEM)  through  a  uniform  electron  cloud 
In  which  no  electron  collisions  occur,  either  with  gas  molecules  or  container 


walls.  The  equation  of  motion  for  the  electrons  gives  the  well-known  dlsperson 

k=k,.  (1-  )  (2-4) 


relation 


I 


where  the  propagation  vectors  Inside  and  outside  the  electron  cloud  are,  respec¬ 
tively,  k=2  jr  f  and  k^*  2  jt  f/p^e^  .  The  plasma  frequency 


f 


N  e‘ 


4  7r  ^  m 


(2-5) 


where  N  Is  the  number  of  electrons  per  m  .  The  electron  cloud  is  in  effect  a 
dielectric  region  with  specific  capacity  K=(k/kQ)  =  e/e^j  .  Use  of  the  appropriate 
boimdary  conditions  for  Incident,  reflected,  and  transmitted  plane  waves  leads  to 
expressions  for  T  and  H in  terms  of  the  plasma  length,  d.  These  are 


r  = 


j/kL 


ik,,d 


(1+/K)*  e^^^^-d-ZiC)^ 


(2-6a) 
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When  the  medium  is  lossless,  we  must  have  I  /^l  in  order  to 


(2-6b) 


conserve 


energy.  If  die  effect  of  elastic  collisions  is  taken  into  account  by  including  the 
usual  Lorentz  damping  term  in  the  equation  of  motion,  the  dispersion  relation  is 
modified  by  the  inclusion  of  a  collision  frequency,  v,  and  equation  (2-4)  becomes 


k=ko  1- 


(2-7) 


Figure  2. 5  is  a  plot  of  |  |  from  (2-6a)  for  several  plasma  lengths  d  expressed  in 

fractions  of  the  plasma  wavelength  (Xp*c/fp).  The  solid  curves  correspond  to 
1/  s  0  and  the  dotted  curves  show  the  effect  of  including  a  particular  collision  fre¬ 
quency.  The  effect  of  the  actual  non-uniformity  in  the  electron  distribution  is 
discussed  Section  2.3.  Since  some  frequency -dependent  fluctuation  will  occur  in 
the  measurements  due  to  line  reflections,  the  minimum  desirable  length  should 
evidently  not  be  less  than  about  The  heater  length  used  for  the  present 

diode  is  31.5cm  and  due  to  the  axial  temperature  profile,  this  restricts  the  active 
region  to  about  23cm.  This  corresponds  to  a  wavelength  of  92cm  or  a  plasma 
frequency  of  326  Me  and  a  density  p  of  1.3  x  10^  electronG^cm^.  Since  this  is 
near  the  density  limit  noted  above,  we  conclude  that  the  span  of  electron  plasma 


frequencies  available  for  good  data  is  restricted  to  the  region  100  to  400  Me. 
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The  composition  and  activation  of  the  cathode  oxide  surface  followed  nor¬ 
mal  procedures  recommended  by  tube  manufacturers.  An  unusual  problem  was 
encountered  in  coating  the  long,  inside  surface.  Ba  CO 2  and  Sr  CO 2  (Sylvania 
CIO  triple  carbonate)  in  a  nitrocellulose  lacquer  and  amyl  acetate  binder  was 
applied  as  a  liquid  coating  on  this  surface  by  the  cathophoretic  process.  A  tech¬ 
nique  of  spin-drying  was  developed  to  promote  a  uniform  depth  of  coating.  The 
coated  cathode  was  then  mounted  in  the  vacuum  system  and  activated  by  heating 
from  the  diode  heater.  Difficulties  in  determining  the  actual  pressure  and  tem¬ 
peratures  at  the  cathode  surface  as  well  as  an  axial  temperature  gradient  were 
encountered;  the  technique  is  still  in  a  formative  stage.  However,  two  cathodes 
have  been  successfully  developed. 

The  Anode  Cooling  System.  The  single-ended  water  cooling  channel  was 
designed  to  operate  at  an  available  tap  pressure  of  50  p.s.i.  For  a  heat  flux  of 
100  watts^cm^  and  the  diode  geometry  already  described,  this  corresponds  to  a 
bulk  temperature  rise  of  30.3*^C  and  a  heat  dissipation  of  3160  watts.  Local 
boiling  temperatures  should  not  yet  be  reached  according  to  conventional  film 
drop  calculation.  Operation  of  electrically  heated  test  models  suggest,  however, 
that  local  boiling  does  commence  near  1500  watts  input  and  the  region  of  instable 
film  boiling  is  reached  above  2500  watts.  Consequently,  an  upper  limit  of  2400 
watts  has  so  far  been  imposed  on  the  diode  operation.  This  corresponds,  by 
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equation  (2-3),  to  an  electron  density  of  1.43x10^  rather  than  the  anticipated 
density  of  1.6x10®.  The  discrepancy  in  heat  transfer  properties  is  not  yet  under¬ 
stood  and,  because  of  Inadequate  test  models,  may  not  actually  exist  in  the  electron 
diode .  The  exceedingly  small  dimensions  for  the  stainless  steel  tubing,  shown  in 
Figure  2. 3,  make  support  brackets  for  the  inner  tube  impractical  and  consequently 
lack  of  concentricity  is  unavoidable.  This  may  provide  an  explanation  for  the 
discrepancy. 

Coaxial  Support  Beads.  The  center  conductor  of  the  coaxial  loop  is  sup¬ 
ported  by  ten  dielectric  beads  each  of  which  presents  a  dielectric  discontinuity 
and  hence  a  possible  reflection.  These  reflections  have  been  effectively  com¬ 
pensated  in  the  low  frequency  region  below  3000  Me  by  adjustment  of  the  bead 
radial  dimensions.  The  six  beads,  shown  in  Figure  2.1,  inside  the  vacuum 
system  are  subjected  to  temperatures  up  to  200°C  .  Consequently,  quartz  (K=4.1) 
was  used  in  this  region.  The  rest  of  the  beads  are  made  of  teflon. 

The  large  teflon  beads  have  a  three -fold  design  requirement  and  were 
selected  after  considerable  e:q)loration.  They  must  provide  good  RF  impedance 
matching,  mechanical  strength  against  accidental  vibrations,  and  a  tension  force 
of  up  to  100  pounds.  They  must  also  insure  a  high  vacuum  seal.  The  beads  were 
surface  etched  to  provide  a  bonding  surface.  A  vacuum -tight  bond  to  the  nickel 
coaxial  walls  was  achieved  by  .010  in.  thick  coating  of  epoxy  resin.  Metal  and 
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teflon  walls  were  scored  0,002  in.  deep  to  promote  bonding  and  the  epoxy  was 
built  up  on  successive  applications  to  seal  against  bubble  leaks.  External  cooling 
coils  insure  that  the  teflon  temperatures  remain  below  lOO^C  .  Operating  e;g>eri- 
ence  with  these  beads  has  been  quite  satisfactory. 

Vacuum  System.  The  vacuum  system  was  designed  to  provide  a  vacuum 
below  10~'^mm  Hg  which  is  sufficient  for  the  protection  of  the  activated  cathode. 
The  vacuum  console  was  engineered  by  Varian  Associates;  most  of  the  components 
are  visible  in  Figure  2.2.  The  150  liter  test  chamber  is  evacuated  throu^  a 
5  in.  manifold  by  two  valved  'Vacsorb"  adsorption  pumps  and  a  "Vacion"  ion 
pump.  A  Phillips -Granville  controlled  leak  valve  also  coimects  to  the  manifold. 
All  test  penetrations  into  the  system  are  made  through  the  removable  cover  plate 
by  means  of  several  removable  port  flanges.  These  penetrations,  all  insulated, 
include  the  two  coaxial  seals,  six  thermocouple  leads,  and  five  heater  and  diode 
current  leads.  The  test  chamber  walls  can  be  preheated  to  90°C  for  outgassing 
and  are  water-cooled  during  operation.  Starting  pressures  are  indicated  by  a 
thermocouple  vacuum  gauge  and  high  vacua  by  the  ion  pump  current  directly.  The 
ion  pump  has  a  flat  pump  rate  characteristics  of  140  liters /sec  for  air;  the  pump 
rate  varies  with  the  kind  of  gas  being  pumped. 

Operating  experience  has  shown  that  when  the  chamber  is  outgassed  and 
without  the  diode  heater  on,  indicated  pressures  of  4x10'^  mm  Hg  are  normally 
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obtained,  and  this  can  be  accomplished  from  open  tank  conditions  in  about  12 
hours.  With  the  heater  on,  the  level  rises  to  as  high  as  4x10"'^ mm  Hg.  These 
pressures  are  believed  to  be  quite  satisfactory  for  a  replaceable  cathode.  Longi¬ 
tudinal  slits  in  the  connecting  tapered  sections  provide  the  pumping  orifice.  Thus 
there  is  some  uncertainty,  due  to  the  closed  nature  of  the  diode  design  of  the 
actual  pressure  in  the  diode  during  activation  outgassing.  This  has  hampered 
the  activation  operations. 

I 

I 

RF  System.  The  microwave  detection  is  accomplished  with  conventional 
apparatus.  Figure  2.  6  is  a  typical  arrangement  used  for  measuring  both  trans¬ 
mitted  and  reflected  power  relative  to  input  power.  If  the  reflections  from  the 
coaxial  line  without  electrons  is  small,  the  measurement  of  the  transmission 
coefficient  is  esjjecially  simple:  | '?'!  is  just  the  ratio  of  transmitted  voltage  ampli¬ 
tudes  for  the  heated  electron  tube  with  plate  voltage  on  relative  to  plate  voltage  off. 
An  alternative  method  for  measuring  reflected  power  utilizes  a  slotted  line  in 
place  of  the  directional  coupler  and  the  VSWR  is  found  with  plate  voltage  on  and  off. 
In  the  reflection  case,  losses  in  the  cable  and  the  tube  must  be  accounted  for. 

The  performance  of  the  line  and  necessary  corrections  to  the  data  are  discussed 
in  Section  2. 3. 

2.  3.  Experimental  Results  and  Analyses 

In  the  schematic  of  the  experiment  (Figure  2. 1),  the  RF  input  and  output 


23 


-  THE  UNIVERSITY  OF  MICHIGAN  - 

4386-2-F 

lines  are  shown  eneterlng  the  vacuum  console  through  the  cover  plate.  The  block 
diagram  of  the  microwave  circuit  is  shown  in  Figure  2.  6,  The  power  reflection 
and  transmission  coefficients  of  the  circuit  between  the  input  and  the  output 
reference  planes  are  shown  in  Figure  2. 7.  This  measurement  is  carried  out  with 
the  cathode  hot  but  with  no  DC  voltage  on  the  diode.  The  reflections  in  the  loop 
are  small  from  100  to  1000  Me,  moderate  from  1000  to  3000  Me,  and  unsatis¬ 
factorily  large  from  3000  to  4000  Me.  The  power  transmission  coefficient  varies 
from  0. 9  at  200  Me  to  0.  6  at  2000  Me;  in  general  it  is  decreasing  with  increasing 
frequency. 

This  variation  is  due  primarily  to  the  wall  losses  of  the  coaxial  line. 

From  the  behavior  of  the  transmission  and  the  reflection  parameters  one  may 
conclude  that  the  microwave  structure  of  Figure  2.1  may  be  represented  by  lossy 
transmission  lines  connecting  the  diode  to  the  input  and  the  output  reference  planes, 
as  shown  in  Figure  2.  6.  This  is  a  sufficiently  good  approximation  for  frequencies 
from  100  to  3000  Me.  For  such  a  circuit  the  power  transmission  coefficient  of 
the  diode  is  readily  obtained  by  monitoring  the  power  at  the  output  terminal  with 
the  cathode  at  the  proper  temperature  and  the  DC  voltage  on  the  diode  alternately 
on  and  off.  The  ratio  of  the  power  output  with  the  DC  voltage  on,  to  that  with 
the  voltage  off,  is  then  the  power  transmission  coefficient  for  that  section  of  the 
coaxial  line  which  contains  the  electrons. 
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FIGURE  2.  6:  TYPICAL  DETECTION  SYSTEM 
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The  measurement  of  the  power  reflection  coefficient  of  the  diode  is  not  as 
simple.  The  power  reflection  coefficient  is  first  measured  at  the  input  terminal, 
and  then  transferred  to  the  input  of  the  coaxial  diode.  Obviously  for  this  operation 
one  must  know  the  losses  in  the  connecting  line.  The  measured  results  of  the 
attenuation  (with  the  plate  voltage  off)  between  the  input  and  the  output  terminals 
(i.  e.  for  the  entire  circuit)  are  shown  in  Figure  2.8.  The  lower  curve  is  com¬ 
puted  from  the  skin  depth  loss  formula  by  assuming  reasonable  wall  conductivities 
and  permeabilities.  It  exhibits  the  expected  dependence  on  the  square  root  of  the 
frequency.  A  curve, which  is  a  reasonable  fit  to  the  measured  points  shows  the 
same  frequency  dependence,  but  the  actual  losses  are  about  1.5  times  higher. 

This  is  not  unusual  in  view  of  some  uncertainty  of  the  wall  conductivity  and  es¬ 
pecially  of  the  permittivity  of  the  nickel  walls.  The  appropriate  fraction  of  this 
loss  was  used  to  advance  the  refelection  coefficient  from  the  input  terminal  to  the 
input  of  the  diode. 

The  nominal  length  of  the  first  cathode  activated  in  the  system  was  14cm. 
The  power  reflection  and  transmission  coefficients  of  this  cathode  as  a  function  of 
frequency  for  anode  currents  of  0.5,  1.0  and  2.0  amperes  (that  were  obtained  at 
nominal  anode  voltages  of  80,  120  and  190  volts  respectively)  are  shown  in  Figures 
2.  9,  2. 10  and  2. 11.  The  nominal  length  of  the  second  cathode  was  20.  3cm.  The 
current  densities  were  selected  to  correspond  as  closely  as  possible  to  the  res¬ 
pective  runs  of  the  first  cathode.  This  data  is  presented  in  Figures  2. 9a,  2. 10a 
and  2. 11a. 
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FIGURE  2.  8:  LOOP  ATTENUATION  LOSS  VS  FREQUENCY 
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FIGURE  2.  9:  DIODE  POWER  REFLECTTION  AND  TRANSMISSION  COEFFICIENTS 
VS.  FREQUENCY  AT  0.  5  AMP  ANODE  CURRENT  (14  cm  CATHODE) 
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FIGURE  2. 10:  DIODE  POWER  REFLECTION  AND  TRANSMISSION  COEFFICIENTS 
VS.  FREQUENCY  AT  1  AMP  ANODE  CURRENT  (14  cm  CATHODE) 
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FIGURE  2. 11:  DIODE  POWER  REFLECTION  AND  TRANSMISSION  COEFFICIENTS 
VS.  FREQUENCY  AT  2  AMPS  ANODE  CURRENT  (14  cm  CATHODE) 
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Note  that  good  repeatability  is  obtained  between  the  three  or  four  runs 
shown.  The  accuracy  of  the  measurements  is  of  the  order  of  1 5  percent.  The 
interesting  features  are,  1)  the  power  transmission  coefficient  is  flat  over  an 
appreciable  frequency  Interval  at  lower  frequencies,  and  2)  there  is  an  absorption 
peak  at  a  high  frequency,  and  3)  the  power  reflection  coefficient  remains  small 
at  all  times.  For  a  lossless  medium  the  power  reflection  transmission  coef¬ 
ficients  should  add  up  to  unity  at  all  frequencies.  Since  this  does  not  happen,  it 
is  clear  that  the  electron  stream  appears  as  a  lossy  medium  to  the  microwaves 
even  up  to  fairly  high  frequencies. 

In  the  coaxial  configuration  the  interaction  of  the  space  charge  limited 
electron  stream  with  the  incident  electromagnetic  field  is  a  reasonably  well 
defined  boundary  value  problem.  From  the  solution  of  Maxwell's  equations  sub¬ 
ject  to  the  appropriate  boundary  conditions  one  can  predict  the  reflection  and  the 
transmission  coefficients  for  the  diode.  However,  the  complete  solution  is  of 
considerable  difficulty  and  it  still  remains  to  be  worked  out.  Only  an  approximate 
solution  will  be  attempted  here.  The  metal  walls  of  the  diode  appear  to  the 
electrons  as  large  molecules,  and  hence  we  may  talk  of  electron-molecule  col¬ 
lisions.  The  collisions  between  the  electrons  and  the  hypothetical  molecules  are 
inelastic  whereby  the  total  momentum  of  the  electron  is  reduced  to  a  low  value. 

DC  field  intensity  is  sufficiently  strong  so  that  the  drift  velocity  of  the  electrons 
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is  large  compared  to  the  random  velocity.  The  cathode  to  anode  spacing  consti¬ 
tutes  the  mean -free -path  distance  and  the  inverse  of  the  DC  electron  transit  time 
may  be  taken  as  the  collision  frequency  (v).  The  collision  frequency  absorption 
mechanism  is  normally  derived  for  elastic  collisions.  Here  we  are  arbitrarily 
extending  the  mechanism  to  a  hypothetical  gas  with  only  inelastic  collisions 
allowed. 

As  an  example  the  plasma  frequency  in  a  15  cm  diode  is  given  in  Figure 
2.12.  These  curves  are  computed  from  the  electron  densities  obtained  from  the 
Child's  law  potential  distribution  in  the  diode.  Using  the  hypothetical  plasma  pic¬ 
ture  we  may  say  that  the  plasma  frequency  variation  radially  in  the  interelectrode 
space  is  a  microscopic  behavior,  because  it  is  on  a  scale  smaller  than  the  electron 
mean-free-path.  In  order  to  use  the  equivalent  transmission  line  analysis  we 
have  to  obtain  a  macroscopic  parameter,  i.e.  an  equivalent  plasma  frequency. 

If  one  approximates  the  plasma  frequency  in  the  diode  by  a  step  function,  the 
lower  value  (f  2)  given  by  the  flat  portion  of  the  curves  and  the  higher  value  (f  1), 
describing  the  conditions  near  the  cathode,  given  by  the  location  of  the  first  absor¬ 
ption  peak  in  the  transmission  curves,  then  one  may  obtain  an  equivalent  plasma 
frequency  from  the  solution  of  the  principal  mode  (Klopfenstein,  1954).  An  alter¬ 
nate  treatment  is  attempted  in  Appendix  A  of  this  report.  Since  the  arguments 
of  the  characteristic  equation  for  the  principal  mode  in  the  region  of  interest  are 


( 
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FIGURE  2. 12:  PLASMA  FREQUENCY  VS.  DIODE  RADIUS  FOR  DIFFERENT 
ANODE  CURRENTS 
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considerably  less  than  unity  we  may  approximate  both  solutions  of  the  Bessel's 
equation  and  of  the  modified  Bessel's  equations  by  the  first  terms  in  the  series. 


This  gdves  the  plasma  frequency  in  a  very  simple  form: 


r  AUr  /r  ) 


if  /f  )*  =  1  - - tils 

P  1+  i-gjtf 


(2-8) 


l-(fi/f)*  JUtirfy/rg) 
where  r^  denotes  the  position  of  the  step  in  the  electron  density,  and  f  is  the 
frequency  of  the  incident  microwave  radiation.  The  values  of  the  parameters 


used  in  the  subsequent  calculations  are  summaried  in  Table  2.1, 


TABLE  2.1  DIODE  PARAMETERS 


Figure 

Anode 

^2 

v/2ir 

Cur  rent(  amps) 

(cps) 

(cps) 

(l/sec) 

2.9 

0.50 

1.72x10® 

12.25x10® 

1.11x10® 

2.  9a 

0.83 

1.79x10® 

12.25x10® 

1.11x10® 

2.10 

1.00 

2.09x10® 

15.0x10® 

1.35x108 

2.10a 

1.53 

2.13x10® 

16.0x10® 

1.35x10® 

2.11 

2.00 

2.66x10® 

19.5x10® 

1.65x10® 

2. 11a 

3.05 

2.70x10® 

21.0x10® 

1.65x10® 

r^=  0. 198cm,  rp=  0. 442cm,  ^®2cm 

Having  determined  the  effective  plasma  frequency  fp,  and  the  collision 
frequency  v  for  the  diode  we  may  proceed  to  use  the  standard  transmission  line 
formulas  presented  in  the  previous  section.  The  voltage  transmission  and 
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reflection  coefflcientSj  respectively,  of  such  a  medium  of  length,  d,  are  given 
by  equations  (2-6a)  and  (2-6b).  The  squared  amplitudes  of  (2-6a)  and  (2-6b)  for 
the  appropriate  parameters  given  in  Table  2,1,  are  shown  either  as  solid  or 
dashed  curves  in  Figures  2.  9  to  2.11a  .  The  computed  power  transmission  coef¬ 
ficient  predicts  absorption  as  the  microwave  frequency  approaches  the  hi^ 
plasma  frequency  close  to  the  cathode.  As  f  decreases  sufficiently  below  f^, 
then  the  right-hand  side  of  (2-8)  becomes  negative  (i.e.  the  effective  plasma 
frequency  for  the  diode  is  imaginary)  and  this  together  with  the  collision  fre¬ 
quency  predicts  a  growing  wave.  Therefore,  the  power  transmission  coefficient 
swings  above  unity.  The  power  reflection  coefficient  in  this  area  shows  a  spike 
approaching  unity.  As  the  microwave  frequency  approaches  the  lower  plasma 
frequency,  f2,  the  power  reflection  coefficient  increases  and  the  power  trans¬ 
mission  coefficient  decreases.  However,  the  changes  are  not  rapid  enough  to 
give  a  good  fit  to  the  measurements.  For  frequencies  considerably  above  the 
lower  plasma  frequency  the  absorption  is  higher  than  predicted  by  the  theory. 

For  these  frequencies  the  transmission  and  reflection  behavior  is  only  slightly 
influenced  by  the  high  plasma  frequency  near  the  cathode.  The  calculated  power 
transmission  curves  level  out  only  below  100  Me. 

The  above  model  predicts  the  behavior  at  the  high  frequencies  reasonably 
well,  including  the  absorption.  If  a  better  approximation  could  be  obtained  which 
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incorporated  the  gradually  changing  character  of  the  plasma  frequency  near  the 
cathode,  probably  one  could  account  fully  for  the  absorption  dip.  The  main  dif¬ 
ficulty  is  at  the  frequencies  below  fi  .  The  predicted  attenuation  effects  do  not 
fall  off  fast  enough  to  match  the  experimental  curves.  There  does  not  appear  to 
be  any  collision  frequency  which  would  predict  this  behavior. 

Since  the  electron  density  variation  in  the  inverted  diode  is  relatively  flat 
and  since  the  first  higher  order  mode  occurs  around  15  Gc,  the  solution  of  the 
transmission  and  reflection  coefficients  on  the  basis  of  the  principal  modes  only 
is  a  reasonably  good  approximation.  If  the  electron-wall  interaction  truly 
simulates  inelastic  gas  collisions,  then  the  data  indicate  that  absorption  from 
such  collisions  is  greater  than  from  the  normal  elastic  collisions  for  a  sub¬ 
stantial  frequency  interval. 
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number  density  of  gas  particles 
gas  pressure 

radial  coordinate  (of  cylindrical  coordinates) 

anode  radius 

cathode  radius 

gas  temperature 

anode  temperature 

cathode  temperature 

ion  drift  velocity 

electron  drift  velocity 

a  -  number  of  ionizing  collisions  suffered  by  one  electron  per  unit 
path  length  in  the  field  direction  (first  Townsend  coefficient) 

number  of  recombinations  per  second  per  unit  volume  per  unit 
density  of  electrons  per  unit  density  of  ions 

permittivity  of  free  space 

fraction  of  the  kinetic  energy  of  an  ion  \^ich  it  loses  on  the  average 
in  an  elastic  collision  with  a  gas  molecule 

fraction  of  the  kinetic  energy  of  an  electron  which  it  loses  on  the 
average  in  an  elastic  collision  with  a  gas  molecule 

ion  mean  free  path  for  elastic  collisions 

electron  mean  free  path  for  elastic  collisions 
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=  positive  charge  density 
p  =  magnitude  of  the  negative  charge  density 

=  total  elastic  collision  cross  section  of  the  g^s  particles  for  ions 
a  =  total  elastic  collision  cross  section  of  the  gas  particles  for  electrons 

Wherever  the  vector  quantities  7+  and  ^  are  written  without  the 
arrow,  their  absolute  value  is  referred  to. 

In  the  first  part  of  this  Section,  the  case  of  no  ionization  is  treated.  Since 
no  distinction  between  electrons  and  ions  has  to  be  made  in  that  part,  the 
quantities  j,m,v,  x,X,p,g,  are  used  without  the  minus  sign  in  the 
subscript.  It  is  understood  that  these  quantities  refer  to  electrons. 
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ni 

THE  ELECTRON  DISTRIBUTION  FUNCTION 

3, 1  Introduction 

In  order  to  deduce  from  experiment  the  correct  relation  between  the  electric 
properties  of  a  plasma  and  its  electron  distribution,  it  is  necessary  to  use  an 
independent  method  for  the  determination  of  the  distribution.  The  derivation  of 
an  expression  for  the  electron  density  in  terms  of  directly  measurable  observables 
is  based  on  the  fact  that 

T  =  -  P  V  ,  (3-1) 

where  "j  =  current  density, 

p  =  electron  charge  density, 

V  =  electron  drift  velocity, 

and  the  further  fact  that  the  drift  velocity  is  determined  by  the  gas  temperature 
and  pressure,  the  electric  field,  and  the  collision  cross  section  of  the  gas  mole¬ 
cules  for  electrons.  The  current  density  can  easily  be  expressed  in  terms  of 
the  total  current  which  is  measurable.  The  electric  field  can  be  determined 
through  Poisson's  equation  after  the  charge  density  has  been  eliminated  from  it, 
and  the  other  quantities  can  all  be  measured  directly  so  that  the  electron  density 

-V 

can  then  be  calculated.  The  minus  sign  in  (3-1)  is  the  result  of  the  fact  that  j  is 
in  the  direction  of  the  conventional  current,  while  "v  is  in  the  direction  of  flow  of 
the  electrons,  p  is  then  defined  as  the  absolute  value  of  the  electron  charge 
density. 
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3.2  Case  of  No  Ionization 

To  begin  with  the  case  of  no  ionization  will  be  considered.  Furthermore, 
it  will  be  assumed  that  the  only  collisions  which  the  electrons  undergo  are  elastic. 
In  that  case  the  electron  drift  velocity  is  given  by  the  expression 

where  e  =  electronic  charge, 
m  =  electronic  mass 

k  =  fraction  of  the  kinetic  energy  of  an  electron  which  it  loses  on 
the  average  in  an  elastic  collision  with  a  gas  molecule, 

E  =  magnitude  of  the  component  of  the  electric  field  in  the  direction  of 
the  electron  drift  (radial  direction  in  the  present  case), 

X  =  electron  mean  free  path  for  elastic  collisions. 

The  mean  free  path  depends  on  the  total  elastic  collision  cross  section  a  of 
the  molecules  for  the  electrons  in  the  following  way 

X 

where  n  is  the  number  density  of 

where  T  is  the  absolute  temperature  of  the  gas  and  k  is  Boltzmann's  constant. 
Substituting  equations  (3-3)  and  (3-4)  into  (3-2)  one  obtains  the  expression 

^  See  for  example  equation  (4. 21)  in  von  Engel  (1955)  . 


no 


(3-3) 


gas  molecules.  But  the  equation  of  state  is 
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In  the  present  investigation  the  temperature  is  a  function  of  position,  since 
a  hot  cathode  discharge  is  being  used.  Let  the  cathode  radius  and  temperature 
be  r  and  T  respectively  and  the  corresponding  values  for  the  anode  r^  and  T.  . 

C  C  cl 

If  the  outer  cylinder  is  the  cathode,  then  r^  >  r^^.  In  the  steady  state  the  tempera¬ 
ture  distribution  inside  the  gas  is  given  by  the  time-independent  Fourier  law. 
which  is  in  cylindrical  coordinates 

where  the  heating  of  the  gas  due  to  electron-molecule  collisions  has  been  ne^ected 
and  the  temperature  has  been  considered  a  function  of  r  only.  This  latter  assum¬ 
ption  imidies  that  the  cylinder  is  taken  to  be  of  infinite  length  since  the  variation 
of  T  with  respect  to  z  near  the  bases  of  the  cylinder  has  been  ne^ected.  The  gas 
temperatures  at  the  radii  r^  and  r^^  are  T^  and  T^  respectively.  Solution  of 
equation  (3-6)  for  those  boundary  conditions  gives  the  following  temperature  dis¬ 
tribution: 

T  =  -^7  <Tclog^  -  T.log  ^  )  ,  (3-7) 

log  - 

where 

If  this  expression  for  the  temperature  is  substituted  into  (3-5),  the  drift  velocity 


4386-2-F 


becomes 


v  =  c/yogfj;;-Talog 


-L  /T 

'c  /  P  ' 


where 


(3-8) 

(3-9) 


a  log 


Is  a  constant. 

Noting  that  the  current  density  and  radius  vectors  are  parallel  we  have 


I  =  2t/r  •  =  2irjf  T  j 


^  2  wrt  2  ir  r  ’ 


(3-10 


where 


I  current  in  external  circuit 
J  =  ill  -  current  per  unit  length  of  cylinder, 
I  -  length  of  cylinder. 


and  J  and  r  are  the  magnitudes  of  j  and  r  respectively.  Scalar  multiplication 
of  (3-1)  by  2  t"?  yields 


2irrj=J  =  2?rrpv,  (3-11 

where  the  minus  sign  on  the  right  hand  side  disi4>pears  because  the  vectors 
and  "y  have  opposite  directions.  Substituting  (3-8)  into  (3-11)  one  obtains 

J _ 

I.  ^ 

.  (3-lJ 


P  =■ 


2t  rc 


i/f 
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Hence  Poisson's  equation  becomes  in  terms  of  the  electric  field 


-  -^{r  E)  = 
r  dr 


2  J 


r.ocJvogi-Talogi 


where  is  the  permittivity  of  free  space  and  the  minus  sign  results  from  the 
fact  that  p  is  the  absolute  value  of  a  negative  charge  density.  Carrying  out  the 
differentiation  and  regrouping  terms  yields 


2  J/y 


Now  let 


=  Y. 

logr  =  z, 


logra  = 
logr^  = 


*a' 


Then  the  above  differential  equation  becomes 


dY^3 
dz  2 


3  J  /p* 


c  Jt  (z-z  )  -  T  (z-z  ) 
o  »  c  a  a  c 


=  0 


with  the  boundary  condition 


Y  =  0  at  z  =  z 
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and  the  minus  sign  on  the  right  hand  side  has  been  absorbed  by  interchanging  the 
limits  of  integration.  Hence,  the  electric  field  as  a  function  of  radius  is 


E 


c 

0 


6J  ^ 

c  /t  -  T' 
'  c  a 


(r  )2 
c 


T 

c 


T  -  T 
c  a 


Ta 


T  -  T 
c  a 


3 

,  (3-13) 


where 


Since  the  integral  appearing  in  (3-13)  cannot  be  evaluated  in  closed  form,  either  a 
numerical  integration  has  to  be  performed,  or  the  integrand  has  to  be  {4>proximated. 
Two  methods  of  approximation,  one  for  large  and  one  for  small  xi  will  be  briefly 
discussed  in  the  appendix  to  this  section. 

We  are  now  in  a  position  to  write  down  the  expression  for  the  charge  density. 
Substituting  equation  (3-13)  into  equation  (3-12)  and  simplifying  the  resulting 
expression,  one  finally  obtains 
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(3-14) 


Equation  (3-14)  expresses  the  charge  density  in  terms  of  either  directly 

measurable  observables  or  fundamental  constants  except  for  the  two  quantities 

X  and  a  ,  which  have  to  be  obtained  phenomenologically  from  published  data. 

X  for  several  gases  may  for  example  be  obtained  from  Fig.  63  of  vonEngel  (1955), 

E 

where  the  constant  values  of  X  for  low  —  are  applicable  here  since  they  refer  to 
the  case  of  no  inelastic  collisions,  which  has  also  been  assumed  in  the  present 
investigation. 

3.3  Introduction  of  Ionization 

In  order  to  take  the  effect  of  ionization  into  account,  the  positive  and 
negative  charge  densities  and  p  are  introduced.  The  corresponding  current 
densities  are  and  j  .  The  positive  (negative)  sign  in  the  subscripts  here  and  in 
other  quantities  implies  that  said  quantities  refer  to  ions  (electrons).  Although 
the  gas  as  a  whole  may  be  neutral,  this  is  not  true  of  the  local  charge  density 
because  of  the  d.c.  field.  The  electric  field,  charge,  and  current  densities  are 
related  by  the  following  equations 
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frE)  =  f-  (P+-P.>.  0-15) 

0 

-^P+P-)  =  ».  0-16) 

-  r  =  0.  (3-17) 

where  "  P+  =  magfnitudes  of  ion  (electron)  current  densities, 

=  number  of  ionizing  collisions  suffered  by  one  electron  per  unit 
path  length  In  the  field  directions  (first  Townsend  coefficient), 

=  number  of  recombinations  per  second  per  unit  volume  per  unit 
density  of  electrons  per  unit  density  of  ions, 

and  V,  are  the  drift  velocities  of  the  ions  and  electrons. 

Equations  (3-16)  and  (3-17)  express  the  conservation  of  electrons  and  ions  respect¬ 
ively.  The  ionization  coefficient  depends  on  the  ionization  cross  section  and  the 

gas  density.  Because  of  the  dependence  of  the  ionization  cross  section  on  electron 

E 

energy,  is  also  a  function  of  electron  energy  and  thus  of  —  .  It  will  be  assumed 
that  recombination  is  negligible  compared  with  ionization.  This  is  reasonable, 
since  ionization  is  causbd  by  electron-neutral  collisions,  while  recombination  is 
generally  caused  by  three-body  electron-neutral-ion  collisions,  and  the  degree  of 
ionization  is  considered  low  in  the  present  treatment.  Hence,  equations  (3-16)  and 
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(3-17)  simplify  to 

r  +j_+ofrj  =0,  (3-16a) 

r  -  o  r  j_  =  0  ,  (3-17a) 

where  now  the  subscript  i  on  or  has  been  dropped. 

The  physical  picture  described  by  equations  (3-16a)  and  (3-17a)  is  that  of  a 
stream  of  electrons  traveling  from  cathode  to  anode  and  causing  an  avalanche  due 
to  ionization.  The  current  at  the  anode  is  entirely  an  electron  current.  The 
current  at  the  cathode  has  both  an  ion  and  an  electron  component.  The  number  of 
electrons  emitted  by  the  cathode  per  second  plus  the  number  of  ions  hitting  the 
cathode  per  second  and  becoming  neutralized  equals  the  number  of  electrons  which 
are  incident  at  the  anode  per  second.  Hence,  the  difference  between  the  electron 
currents  at  the  anode  and  at  the  cathode  equals  the  ion  current  at  the  cathode, 
which  in  turn  equals  the  number  of  ionizations  per  second  in  the  entire  gas,  since 
ions  do  not  accumulate  and  diffusion  to  the  cathode  is  assumed  to  be  the  only 
mechanism  for  the  removal  of  ions. 

The  solutions  of  equations  (3-16a)  and  (3-17a)  for  the  boundary  conditions 


and 


2  r  r  ^  j  (r  )  =  I 
a  -  a 

4  “ 
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are 


j  -a(r-r  ) 


=  w  [*- 


-air  -  r  )] 

ft 


(3-18) 


(3-19) 


where  I  is  again  the  current  in  the  external  circuit. 

In  the  integrations  of  equations  (3-16a)  and  (3-17a)  a  has  been  considered  constant. 
This  is  not  strictly  true  since  the  electric  field  and  therefore  also  the  electron 
energy  are  functions  of  position  and  or  is  a  function  of  electron  energy.  For  the 
elementary  treatment  presented  here  a  suitably  chosen  average  value  of  a  should 
however  be  satisfactory.  Equation  (3-18)  describes  an  electron  avalanche.  To 
see  this,  one  merely  has  to  write  down  the  ratio  of  the  electron  current  at  the 
anode  to  that  at  the  cathode,  which  is 

I-(ra)  2x  rafjjrg)  ""^*’a"*’c^ 

IJr^j)  2vrj.iJ_(r^.) 

Hence, 

where  d  is  the  radial  distance  from  the  cathode  to  anode. 

It  will  now  be  assumed  that  the  electron  drift  velocity  is  determined  solely 
by  the  electric  field  and  the  elastic  electron-neutral  collisions.  This  is  reasonable 
since  the  average  electron  energy  is  considerably  below  the  ionization  and  even 
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excitation,  potentials.  Therefore  only  the  very  small  fraction  of  electrons  in  the 
high  energy  tail  of  the  distribution  will  undergo  inelastic  collisions,  and  the 
overwhelming  majority  of  electrons  will  suffer  elastic  collisions  only  .  The  same 
is  of  course  true  of  the  ions.  The  drift  velocities  of  ions  and  electrons  are  there¬ 
fore  given  by  equation  (3-5),  viz 


(3-20) 

(3-21) 


where  m^  is  the  mass  of  the  ions  and  X  and  X  have  been  defined  following  equation 
(3-2).  Substituting  the  above  relations  for  the  drift  velocities  into  the  expressions 
for  the  current  densities,  one  obtains 


(3-22) 


where  (3-3),  (3-4)  and  (3-7)  have  been  used  and 


(3-23) 


are  the  total  elastic  collision  cross  sections  of  the  neutral  particles  for  ions 
and  electrons  respectively.  The  positive  and  negative  charge  densities  follow 
from  equations  (3-18),  (3-19)  and  (3-22)  and  are  given  by  the  expressions 
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J  Ll-e 


-Q<r-ro) 


2  '  r  c^/r^log  Fa-VogF^  ff 


(3-24) 


j  -O^r-Ta) 

J  e  “ 


2jrrc /T^logf^-T^logf^  . 


(3-25) 


where  again 


J  =  l/i  . 


Substitution  of  (3-24)  and  (3->25)  into  (3-15)  3deld8  finally 


(3-26) 


Again  introducing  the  transformation 

E^  =  Y 

one  obtains  a  differential  equation  of  the  form 

^  3  Ly  =  1  bi-tb, 

dr  2  r  2  rr— ; — ■  •» 

r ,  b3logr+b4 


where 


2J^ 


2j(Er 


t-O- 


Tc-Ta  , 


Talogrc-T  logr  . 


(3-27) 


The  general  solution  of  (3-27)  for  the  boundary  condition  Y(r„)=0  is 
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Y  =  E 


^  ^  b  i4b2  e 

^  J 


-cKp-r-) 


bglog  ' 


(F  dp 


(3-28) 


A  remark  concerning  the  sign  of  Y  is  in  order  at  this  point.  Since  p^. 
and  p_  are  the  absolute  values  of  the  charge  densities,  it  is  clear  that  the  positive 
square  root  of  E  enters  into  equations  (3-24)  and  (3-25),  But  a  cursory  inspection 
of  the  integral  appearing  in  (3-28)  will  reveal  that  it  may  be  either  positive  or 
negative  depending  on  the  upper  limit  r,  and  if  the  integral  is  negative,  then 

In  order  to  insure  that  one  obtains  the  positive  square  root  of  E,  one  has 
to  take  the  absolute  value  of  the  integral  when  the  expression  for  >fE^  is  substituted 
into  the  expressions  for  p+  and  p_  . 

The  Integral  in  (3-28)  can  be  written  as  the  sum  of  two  integrals.  The  first 

r  ^ 

of  these  can  again  be  reduced  to  an  integral  of  the  type  J  e  dx  by  the  substitu¬ 
tions  log  p  =  z  and  =  x^  .  The  second  integral  can  frequently  also  be 

reduced  to  such  an  Integral  if  a  power  series  expansion  of  the  exponential  function 
is  practical.  To  show  this,  consider  the  integral 

/P' 


^bglogp+b^* 


;  dp  , 


where  p°  (n  integer)  is  proportional  to  a  typical  term  in  the  expansion  of 
-a(p-ra) 


.  Let 


log  p  =  z  . 


I 
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Then 


Again  let 


b  z+b 
3  4  9 

- 7 - = 


Then 


2n+3.,  2  V  k  \  2n-t-3  °4  f 

-5—  (x^-  b4^3)  -  -7-  153  I 

1  dx  =  2  e  J  < 


2n+3  ^2 


dx  . 


In  many  practical  applications  only  a  few  terms  of  the  exponential  series  are 
required  to  achieve  convergence.  For  example  in  our  investigation  E/p  is  of 
the  order  of  200  volts^m  x  mm  From  Figures  4. 34a  and  4. 35  of  Brown 
(1959)  it  follows  that  the  corresponding  value  of  a/p  for  both  helium  and  nitrogen 
is  of  the  order  of  l(cm  x  mm  Hg)  .  Since  a  pressure  of  1mm  Hg  has  been  used, 
a  =  1  cm~^  in  both  cases.  Since,  furthermore  p  -r  varies  from  0  to  0. 245cm, 

cl 

the  e^qxment  a^p-r^)  is  obviously  siifficiently  small  for  convergence  to  occur 
after  only  a  few  terms  depending  on  the  accuracy  desired. 

The  electron  charge  density  is  again  obtained  by  substituting  the  expression 


for  E  from  (3-28)  into  (3-25)  .  After  simplifying  this  yields 

1^ 

^e 


I 

-T 


.A 

’-"2x  [3/V 


-a(r-ra) 


y'<T„io6i^-yog|„) 


c  -(c_^+c 

— yPdp 


(3-29) 
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Similarly  the  positive  charge  density  is 
c  I* pc 


V2jr 


3i2c2 


1  -  e 

yr<T^logf^-T^logfj 


c  -(c_j^+c 

/yog|^-T__iog|^ 


/pdp 


,, 

(3-30) 


It  is  easy  to  check  that,  if  a=0  (no  ionization),  then  p^=0  and  p  goes  over  into 
the  relation  (3-14)  after  the  appropriate  change  in  the  variable  of  integration. 

The  basic  jdiysical  assumptions  made  in  the  derivation  of  the  above  expres¬ 


sions  for  the  charge  densities  are  summarized  below. 


1)  Electron-ion  recominations  are  negligible 

2)  The  ionization  coefficient  a  used  in  this  treatment  is  a  suitably  chosen 
mean  value. 

3)  The  electron  drift  velocity  is  determined  only  by  the  balance  between 
the  elastic  electron-neutral  collisions  and  the  acceleration  of  the  electrons  by 
the  electric  field.  Since  the  drift  velocity  is  the  average  of  the  electron  velocity 
components  in  the  direction  of  the  electric  field,  this  assumption  means  that  in 
computing  this  average  the  distribution  is  cut  off  at  the  minimum  energy  necessary 
to  cause  excitation.  The  assumption  is  therefore  valid  only  if  the  mean  energy  of 
the  electrons  is  considerably  below  the  lowest  excitation  potential  of  the  atoms 

(or  molecules)  in  question  since  in  that  case  only  the  high  energy  tail  of  the 
distribution  function  would  be  ne^ected. 

4)  The  fractional  change  in  electric  field  over  one  mean  free  path  is  small. 


i.  e. 


k  AA 

E  dr 
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This  is  implied  in  the  derivation  of  the  expression  (3-2)  for  the  drift  velocity. 

5)  The  fractional  energy  loss  per  collision  H  is  not  a  function  of  the  electric 
field.  As  long  as  elastic  collisions  only  are  considered,  this  assumption  is  cor¬ 
rect  as  can  be  seen  from  Figure  63  of  von  Engel  (1955).  In  the  case  of  Inelastic 
collisions  this  assumption  no  longer  holds  true.  But  since  according  to  assum¬ 
ption  3)  the  expression  for  the  drift  velocity  used  here  is  based  on  the  absence  of 
inelastic  collisions,  this  last  assumption  is  therefore  consistent  with  the  rest  of 
the  treatment. 

In  conclusion  a  short  discussion  of  the  expressions  for  the  positive  and 
negative  charge  densities  together  with  a  few  suggestions  for  further  work  will 
be  given. 

First  of  all  it  is  immediately  obvious  that  both  the  expression  for  and 

that  for  p  diverge  as  r->-r  .  This  is  a  consequence  of  the  relations  for  the  drift 
**  c 

velocities  (3-20)  and  (3-21).  According  to  these,  v+  =  0  if  E  =  0,  and  the  electric 
field  has  been  taken  as  zero  at  the  cathode.  Since  neither  of  the  current  densities 
vanishes  at  the  cathode,  the  charge  densities  obviously  become  infinite,  if  the 
drift  velocities  vanish.  The  electrons  naturally  do  not  leave  the  cathode  with  zero 
velocity  but  rather  with  a  Maxwellian  velocity  distribution  which  is  determined  by 
the  cathode  temperature.  Furthermore  the  assumption  of  a  slowly  varying 
electric  field  is  especially  bad  in  the  region  near  the  cathode.  In  particular  the 
electric  field  does  not  vanish  at  the  cathode  but,  because  of  the  space  charge,  a 
certain  distance  away  from  the  cathode.  The  region  where  the  electric  field 
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vanishes  is  referred  to  as  the  virtual  cathode.  Because  of  the  divergence  of  the 
expressions  for  the  charge  densities  at  the  cathode,  these  expressions  hold 
throughout  the  discharge  (provided  the  stated  assumptions  are  satisfied)  but 
excluding  the  small  region  immediately  adjacent  to  the  cathode.  Additional  work 
is  required  to  generalize  the  expression  for  the  current  density  by  removing  the 
above  mentioned  objections.  A  correct  treatment  would  of  course  result  in  a 
finite  value  for  the  charge  densities  at  the  cathode. 

According  to  the  fourth  assumption  above  the  fractional  change  in  electric 
field  over  one  mean  free  path  must  be  small.  Because  of  the  cylindrical  geometry 
considered  in  the  present  treatment  E  is  obviously  never  constant.  In  a  vacuum 
and  in  the  absence  of  any  discharge  (i.  e.  no  space  charge)  the  electric  field  is 
proportional  to  l/r  and  therefore 


X  ^ 
E  dr 


X 

r 


The  necessary  and  sufficient  condition  for  the  fourth  assumption  to  be  satisfied 
is  therefore  X4<r 

a 

The  total  charge  in  the  gas  is  given  by 

Q  =  Jip^-pJdY, 


(3-31) 


where  the  integration  has  to  be  carried  out  over  the  volume  occupied  by  the 
discharge.  From  (3-29)  and  (3-30)  one  sees  that  the  integrand  (p^-p  )  equals  a 
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but  had  to  be  abandoned  because  of  the  too  great  comidexlty  of  the  resulting 
equations.  However,  it  would  be  possible  to  lump  the  inelastic  effects  into  the 
fractional  energy  loss  X  and  then  to  use  for  x  a  function  fitted  to  the  experimental 
curves  which  are  given  in  Figure  63  of  von  Engel  (1955).  This  function  can  then 
be  substituted  into  the  expression  (3-2)  for  the  drift  velocity.  If  both  and  X 
can  be  represented  by  the  same  function  proportional  to  some  power  of  E/p,  say 
a  (E/p)’”,  where  a  is  a  constant  and  m  is  any  positive  real  number,  then  the 
differential  equation  (3-26)  can  still  be  solved  by  quadrature  after  the  transfor¬ 
mation  Y  =  E*”  haa  been  carried  out.  The  assumption  of  the  same  functional 
dependence  of  and  X  on  E/p  is  however  questionable.  For  while  the  absolute 
energy  loss  suffered  by  an  Inelastically  colliding  electron  is  the  same  as  that 
suffered  by  an  inelastically  colliding  ion  (namely  the  energy  difference  between 
the  final  and  initial  states  of  the  excited  atom),  the  average  energy  of  the  electrons 
is  higher  for  a  given  electric  field  than  that  of  the  ions.  Hence,  the  average 
fractional  energy  loss  of  electrons  for  Inelastic  collisions  is  smaller  than  that 
of  ions.  It  might  be  conceivable  that  the  different  magnitudes  of  and  x  can 

be  taken  into  account  by  introducing  two  different  coefficients  a^  and  a  with 

m 

a^>a  ,  while  taking  the  same  power  m  in  (E/p)  for  both  electrons  and  ions. 

This  would  still  leave  (3-26)  solvable  in  an  elementary  way. 

However,  one  should  never  lose  sight  of  the  basic  limitation  of  the  entire 
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treatment  \v4iich  becomes  especially  severe  for  higher  values  of  E/p.  The  e}q)res- 
sion  for  the  drift  velocity  used  here  is  based  on  sm  "average  electron  velocity"  and 
therefore  completely  neglects  the  electron  velocity  distribution.  The  above  men¬ 
tioned  phenomenological  introduction  of  a  functional  dependence  of  >^on  E/p  should 
be  clearly  understood  with  this  in  mind.  It  is  actually  merely  an  elementary  way 
of  avoiding  the  real  difficulties.  A  consistent  treatment  would  of  course  have  to 
start  out  from  the  Boltzmann  equation. 
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LIST  OF  SYMBOLS  FOR  SECTION  IV 


effective  absorptioa  coefficient 

(cj)  absorption  coefficient  neglecting  stimulated  emission 

K  photon  energy^  =  hw  =  hi/ 

r  classical  electron  radius,  =6^/mc^ 

o 


y  hc/x  Debye 

K^Ki,  K2  modified  Hankel  functions 

p3  p2 

Pi,  Pj  momenta  before  and  after  collision,  multiplied  by  c;  -  K 

El,  E3  energies  corresponding  to  Pi  and  P2;  Ei  -  Eg  =  K 

/u  me* 

a  fine  structure  constant,  =  6*/hc 

^  polarizability  of  the  neutral  atom 


n 

e 


densities  of  electrons,  ions,  atoms,  scatterers 


electronic  chaise,  to  avoid  confusion  with  Naperian  base  e. 
plasma  angular  frequency 
collision  frequency 

certain  sums  over  atomic  quantum  numbers,  defined  in 
(Barasch  et  al,  1962) 
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IV 

FREE- FREE  ABSORPTION  COEFFICIENT  IN  AN  IONIZED  GAS 

4.1  Introduction 

At  frequencies  exceeding  the  plasma  frequency,  the  most  significant  absorp¬ 
tion  process  for  electromagnetic  radiation  in  an  ionized  gas  is  free-free  absorption. 
In  this  section  we  derive  the  attenuation  coefficient  corresponding  to  free-free  tran¬ 
sitions  of  the  electrons  both  in  the  field  of  singly-charged  ions  and  of  neutral  atoms 
or  molecules.  When  the  electron  density,  neutral  densily,  and  chemical  species 
are  known,  these  coefficients  can  be  evaluated  numerically  for  a  Maxwellian  electron 
distribution  of  given  temperature.  They  are  then  found  to  agree  with  the  attenuation 
predicted  by  the  classical  Lorentz  theory  (i.  e.  the  limit  for  zero  magnetic  field  of 


the  Appleton- Hartree  relations.)  This  agreement  is  to  be  expected;  we  shall  discuss 
the  point.  We  also  derive  the  form  of  the  free-free  absorption  coefficient  for  non- 
Maxwelllan  electron  distributions,  in  terms  of  the  distribution  function. 

4.2  Principle  of  Detailed  Balance  I  Its  Utility 


The  principle  of  detailed  balance  states  that  for  a  physical  system  in  equil¬ 
ibrium,  such  as  a  partially  ionized  gas,  if  we  focus  our  attention  on  a  pair  of  states, 
(1  and  2)  of  this  system,  we  will  find  the  rate  (number  per  second)  of  transitions 


l'^2,  Tj^,  equal  to  the  transition  rate  for  2-^1,  T^j.  (This  is  stronger  than  saying 


^  T  ^) .  Thus,  since  cross  section  for  atomic  processes  are  functions 

J 


66 


THE  UNIVERSITY  OF  MICHIGAN 

4386-2-F 


of  atomic  constants  rather  than  of  the  gas  parameters  such  as  temperature,  we  may 
take  a  situation  in  which  we  know  just  what  processes  are  responsible  for  the  tran¬ 
sitions  in  question,  and  take  a  conceptual  equilibrium  situation.  In  this  equilibrium 
situation,  we  can  write  the  transition  rates  as  a  function  of  the  cross  sections  for 
the  processes  and  their  inverses,  and  by  equating  the  rates,  obtain  a  relation 
between  the  cross  sections  for  the  processes  and  for  their  inverses .  Since  the 
cross  section  is  not  dependent  on  the  state  of  the  gas,  this  relation  is  true  out  of 
equilibrium  too.  Since  in  most  cases  the  cross  section  for  the  inverse  process 
could  be  obtained  directly  only  by  a  very  laborious  quantum-mechanical  calculation, 
the  use  of  detailed  balance  affords  a  great  saving  in  effort. 

Finally,  such  things  as  absorption  coefficients,  which  are  a  function  of  the 
state  of  the  gas,  can  then  be  obtained  as  integrals  of  these  cross  sections  over  the 
distributions  actually  existing. 

In  the  present  case,  we  are  dealing  with  transitions  between  free  states  of 
the  electron  which  are  separated  by  microwave  frequencies,  (free-bound  and 
bound-bound  transitions  are  outside  the  frequency  range  of  interest)  and  we  know 
the  only  operative  processes  are  Bremsstrahlung  (spontaneous  emission),  stim¬ 
ulated  emission,  and  free-free  absorption.  We  shall  therefore  be  able  to  use 
detailed  balance  to  relate  the  cross  sections  for  these  three  processes,  and  then 
formulate  the  effective  absorption  coefficient  in  terms  of  the  previously  obtained 
Bremsstrahlung  cross  section. 
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4.3  Some  Definitions  of  Cross  Sections  and  Coefficients  for  the  Relevant  ProoeBses; 
Their  Relatlms  to  Transition  Rates 

(a)  Bremsstrahluag  1 

The  differential  cross  section  for  Bremsstrahlung  da{K,  P^)  (differential  in/ 
I^oton  energy  K)  describes  emission  of  a  photon  of  enei^  K  when  the  electron 
originally  has  m  Vi  c  =  and  ends  up  with  Pg  such  that  =  2  /iK.  Let  our  gas 

be  composed  of  electrons  and  massive  particles  (ions,  atoms,  or  molecules),  the 
latter  viewed  as  existing  separately,  so  we  consider  cross  sections  for  electron-ion, 
electron-atom,  or  electron- molecule  interactions  with  radiation  separately.  We 
talk  of  states  of  the  system  (electron  +  massive  particle)  and  the  radiation  field. 
Consider  a  volume  of  space  V  for  the  whole  gas.  There  are 

\\ 

V  1  dPi  n.j,  systems/ unit  volume 
Pi 

in  it  with  the  electron  energy  in  dP^  at  P^,  where  n.^  is  the  density  of  massive 
particles.  Since  we  take  the  massive  particles  as  stationary  and  measure  energy 
from  the  zero  of  ionization  potential,  the  energy  of  the  system  is  the  energy  of  the 
free  electron.  Then,  at  last  we  may  say  that  the  number  of  Bremsstrahlung  tran¬ 
sitions  per  second  in  a  volume  of  gas  V  from  the  continuum  range  of  states  with  the 
electron  in  dPi  at  Pi  and  no  photon  to  the  range  of  states  with  the  electron  in  dP2  at 
P2  and  a  photon  with  energy  K  in  dK  (where  -  Pj  =  2  /.(Kand  we  always  consider 
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^  _  energy  loag  In  V/  sec 

energy/photon 


A  dx  I(K)  dK  dn\  „  v 


dn  \ 


ABS' 


(4-4) 


^2 

Finally,  we  put  I  (K)  dK  in  terms  of  the  volume  density  of  photons  n„„  _  (K)dK  = 

I'llUT 

IftOdK 
Ko 


,  and  replace  the  system  distribution  by  the  electron  distribution  multiplied 


by  since  the  massive  particles  are  stationary.  Thus,  we  write  to  define 
®ABS(K,  Pj), 


dn 

Vp.  "  “PHar®-*  "abs  «'  V  • 


(4-6) 


We  also  note  for  later  use  that  the  differential  absorption  coefficient  (differential  in 
electron  energy),  is 


d  Oabs  00=  <to,(P.>o^ABS  «•'■»>• 


(4-7) 


(c)  stimulated  Emission 

This  process  goes  from  the  ranges  dPx  at  P^  down  to  dP2  at  P|,  but  is 
governed  by  the  same  law  as  free-free  absorption.  We  therefore  write  down 
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?  \ 
1 


npHOT<K>dK  C  —  dP^  (K,P,)V  V 


8jrK*dK 


SffPo^dP, 


(4-11) 


“phot®“°15 


«3  - 


dP, 

<hc)^ 


SfK^dK  r  St  d  p7 

V  ^  % 

(he)  (he)"" 


(4-12) 


We  see  now  that  there  is  no  way  of  avoiding  the  usual  choice  n^V  =  1,  i.  e. , 
we  must  work  in  an  "atomic  volume".  Such  a  volume  has  in  it  on  the  average 
n  n_  ^  systems,  which  is  <  1  for  us  when  the  massive  particles  are  neutrals  in  a 

6  X 

weakly  ionized  gas,  so  that  the  situation  is  unsatisfactory,  fluctuations  about  the 
average  being  then  relatively  important.  This  is  the  same  situation  as  met  in  the 
derivation  of  the  Saha  equation,  and,  as  is  customary,  we  adopt  the  procedure 
without  further  comment. 

4. 5  Relation  between  and  Og 


Although  it  does  not  evaluate  these  cross  sections  explicitly,  quantum 
mechanics  indicates  that  they  are  related.  For  the  a  priori  transition  rates  per 


system  between  a  pair  of  states  (that  is,  removing  the  weighting  by  occupation 
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probability)  are  equal  for  free-free  absorption  and  stimulated  emission.  (See 


Schlff  ,  1949),  We  therefore  take  (4-11)  and  (4-12),  remove  from  each  the  statisti- 

— I^/2iikT  ^^e\ 

cal  weighting  factor  e  ^  Implicit  in  the  equate  them.  The 

resulting  relation  is 


or 


4.6  Relatione  in  Equilibrium 

The  statement  of  detailed  balance  is 


Ri  +  Eg  ~ 


Using(4-14)  and  the  fact  that  in  equilibrium, 

8jr  K*  dK 


(K)  dK 

PHOT 


(hc)^ 


while 


dn 

_ S 

dP 


dP.  = 


n  4r  P?  e 
e _ i 


-  I^/2AikT 


(2k  m  kT) 


3/2 


we  obtain 
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(4-14) 


(4-15) 


(4-16) 


(4-17) 
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«oo  -I^/2/ikT 

4t  e 

“  ®  '  (2irn\iTr^ 


.00 


-  n 


dPi 


At  Pi  Pj  e 


-  I^/2jukT 


(4-21) 


/2Tjr 


(2*  fjkT)^^^ 


®ABS  ' 


where  we  have  employed  (4-14)  in  the  second  integral.  But  since  K  is  fixed,  we 
find  from  conservation  of  energy  (which  is  +  2  K  /i)  that 


Pi  dPi  —  Pj  dP2 

We  use  this  to  transform  the  second  integration  to  one  over  P2,  and  obtain 
.00 


(4-22) 


At  P  2 

“if®  =  “e  \  ^,3/2  “abS®  P**' 


-I^/2*ikT 


(4-23) 

Now,  of  course,  is  the  absorption  coefficient  which  would  have  been  obtained  by 
neglecting  stimulated  emission,  or 

-I^/2/4kT 


-00 


4x1^ 


o..  (K)  =  n  \  dPa 
ff  r 


5  e 


(2x  MkT) 


3/2  ABS 


(4-24) 
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We  have  verified^  then,  by  combining  (4-23)  and  (4-24),  that  in  equilibrium 


aJj(K)  =  (l-e"^^*'’^)a^(K)  . 


(4-25) 


In  equilibrium,  the  eiqpliclt  relation  given  by  detailed  balance  between  a*  (K)  and 

dt^Br 

”  (K,  Pi)  ,  obtained  by  combining  (4-18)  and  (4-23),  with  a  change  to  Integrating 

I  ulv 


over  Pi  is. 


,,  .3  .  K/kT 
n  n  (ho)  (e  '  -i 

o'  (K)  = 

2KV  (2v  /ikT)  ' 


(4-26) 


We  use  this  when  (K,  Pi)  is  possible  to  formulate  by  a  direct  procedure,  but 

div  - 

o'  would  be  tedious  to  obtain  by  finding  o .  and  a  directly.  If  (K,  Pi) 

II  Ai>b  OK 

is  g^ven  as  a  function  of  both  Pi  and  P2,  Pj  is  to  be  eliminated  by  expressing  it  as 
y  P?-2K#4'  . 


4.7  Absorption  Coeffloient  out  of  Equilibrium 


The  relations  (4-18)  and  (4-19),  which  are  valid  in  general,  are  the  starting 
points,  to  be  Inserted  into  (4-20).  The  result  is 
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(K,  ^ ) 


8t  K*/i 


da. 


dn^  (Pa,  F)  ^ 


Br 


dK 


(K,  Pi) 


da, 


dn  (Pi,  r) 

G 


BR 


dK 


(K,  Pi) 


/SKm 


(4-27) 


wherein  conserv'  tlon  of  energy  Is  to  be  used  to  eliminate  Pi  in  the  first  integral. 
The  use  of  the  symbol  ?  indicates  the  possibility  of  a  spatially  inhomogeneous 
electron  distribution. 

It  must  be  admitted  that  this  equation  is  of  limited  utility,  since  in  a  non- 

equilibrium  situation,  the  form  of  dn  may  depend  on  the  spectral  density  of 

© 

radiation  present.  However,  in  many  practical  situations,  applied  potential  grad¬ 
ients  or  temperature  gradients  rather  than  externally  generated  electromagnetic 
energy  are  the  essential  factors  in  determining  dn^,  which  may  then  be  obtained 
separately  and  inserted  into  (4-27)  to  give  a  meaningful  result. 

4.8  Summary  of  Expressions  for  the  Bremsstrahlung  Cross-Section 

The  effective  attenuation  coefficient  due  to  free-free  transitions  may  be 

do* 

BR 

evaluated  if  the  Bremsstrahlung  cross-section  —  —  (K,  P^)  is  known.  In  work 

cllv 

under  other  contracts,  expressions  for  this  quantity  have  been  obtained,  which 
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pertain  to  Bremsstrahlung  in  the  fields  of  ions,  atoms,  and  molecules  respectively. 
We  shall  now  list  them. 

In  the  field  of  singly-charged  ions,  we  use  (Barasoh,  I960)  separate  ex¬ 
pressions  for  the  faster  and  slower  electrons.  For  the  faster  ones. 


BR. 
dK  ' 


(K,Pl)-  g  [  2  ^  (Pi  -  Pj)*  +  V* 


2  *7  Pi  Po  I 

[(Pi  +  pj*  +  [(Pi  -  Pj)»  +  7^]  J' 

(4-29) 


while  for  the  slower  ones. 


(4-30) 


where  the  symbol  tT*  denotes  the  first  square  root  written  e^licitly. 

The  criterion  for  distinguishing  between  fast  and  slow  electrons  is  essentially 

that  Pi  >  Pj^,  where  and  where  =  Pj  +  2K/i.  We 

often  find  a  choice  of  P  »  to  be  satisfactory  and  to  simplify  the  subsequent 

1C 

I 

Integrations. 

In  the  field  of  neutral  atoms,  the  Bremsstrahlung  is  treated  by  the  method 
of  Heitler,  using  a  Hamiltonian  which  includes  all  Coulomb  interactions  and  the 
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radiation  field.  The  result  (Barasch  et  al,  1962)  contains  two  contributions,  the 
first  of  which  arises  from  the  static  potential  of  the  atom  and  the  second  from  the 
polarizability.  The  second  should  be  negligible  for  situations  of  practical  interest, 
so  we  write  them  separately  and  n^lect  the  interference  between  them  which  is 
indicated  by  the  formalism.  The  static  contribution  is 


do. 


BR 


dK 


(K,  El)  =  ^  (2  -  K/Ei)  yi  -  K/Ei‘ 


Km 


(4-31) 


Here  a(o)  is  the  elastic  scattering  cross-section  of  the  neutral  atom,  for  low-energy 
electrons  and  Ei  is  the  energy  of  the  Incident  electrons,  or  Ei  -  /2m  . 

The  contribution  due  to  polarizability  of  the  neutral  atom  is  (Barasch  et  al, 

1962) 


(nc; 


1+  a-K/Ei 


1-  /l  -  K/Ei’ 


(4-32) 

Finally,  for  neutral  molecules,  an  expression  identical  to  (4-31)  was  obtained, 
except  that  o(o)  is  now  the  scattering  cross-section  of  the  molecule  for  slow  elec¬ 
trons.  Thus,  again. 


do. 


P-  (K,  El)  =  (2  -  K/Ei)  v/1  -  k/^ 


dK 


(4-33) 
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4. 9  Corresponding  Expresslone  for  the  Effective  Free-free  Absorption  Coefficient 


We  shall  give  the  explicit  forms  assumed  in  equilibrium  by  the  effective 
absorption  coefficient  corresponding  to  these  Bremsstrahlung  cross-sections.  The 
starting  point  Is  equation  (4-26), 

3  K/kT  m 

(K)  \  dP 


-I^/2MkT  „  , 

dPi  Pi  e  ^  (K,  Pi) 


Consider  first  the  free- free  absorption  in  the  field  of  ions.  The  integration 

over  the  lower-velocity  part  of  the  electron  distribution  must  be  performed 

numerically,  and  that  over  the  faster  ones  depends,  of  course,  on  the  value  of  Pi 

chosen  as  the  division  between  slow  and  fast.  If  we  choose,  for  example,  the  lower 

limit  of  integration  as  Pi  -  /jukT ',  and  find  that  the  approximations  Pi  +  P2^^Ti> 

2Kiu 

Pi  Pa  are  valid,  and  also  use  the  relation  Pi  -  P*  =  -  X'  ,  we  find  for  the 

1*1 +  ±'2 

contribution  of  the  faster  electrons 


y  exp  El  (-1/ 2-1/2  ) 


(4-34) 
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For  the  static  contribution  to  free-free  absorption  in  the  field  of  neutral 
atoms,  insertion  of  equation  (4-31)  into  equation  (4-26)  and  changing  the  variable 
of  integration  to  Eg  yields 


a'  (u) 


32  fir  n^n^ocKo) 
3u^  (mkT)^^^ 


kT  J 


dE  E^/^  (E  +  Iko)^/^  (E  +-^) 


(4-35) 


where  we  have  already  taken  h  w/'kT«l,  so  that  1  -  e  h  w/kT. 

We  employ  a  change  of  variable  to  y  =  E/hu,  split  the  integral  into  two 

hu 

terms  resulting  from  the  factor  E  +  ,  and  obtain  integrals  which  are  evaluated 

in  terms  of  the  function  (x)  and  its  derivative.  Recurrence  relations  on  this 
function  are  then  used  to  convert  this  result  to  aKg  (x).  That  is. 


32 


a’  (u)  = 


/F  n^n  o  0(0)  *  ,  , 

'  ae _ hu  .3 


„  3  .  ,  ,,,3/2  kT 
3  w  (mkT)  ' 


[[ 


y  /y  (y+i)  e  + 


+  T  \  dy  yV  (y+i^  e”®^ 


(4-36) 


where  s  =  1kj/2kT  .  Then 
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a'  (w)  »  A  (u,T)  iq  (s/2)^--^  e“/^  Kj  (8/2))  j 

*  A  (w,  T)  e*  iq  (t)  -  e*  iq  (t^j  ,  where  t  =  s/2 

tt  t  t  ■» 

^  iq(t)  +  ^  Ki(t)  ♦  ^  $C^(t)  +  Ka(t))j 
=  A(“<T)  ^  +T(K„  +  K.)]  4- 


A(u,T)  Ka(t) 


Now  for  tiu/2kT  <<1,  Kg  "►  8  »  ■<>  finally  we  may  write 


a*  (u)  “Sf 
static  '  ' 


64  y?  n  a  acr(o)1k  (kT)* 
a  6 


For  the  contribution  of  the  induced  dijxde  In  the  neutral  atom,  we  insert  equation 
(4-32)  into  equation  (4-26).  Thus, 
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16^  It 

a  *  ((j)s  iJk  ^ 

dipole  3,!aPkT 


.  A.  .  (4-39) 

.1- /!-!((£,■ 

Rationalizing  the  denominator  of  the  argument  of  the  logarithm,  we  may  write  the 


integral  as 


^  =  2  J  dE.e-=*/«ip[/f 


(4-40) 


We  now  Introduce  the  substitution 
X  =  “  cosh”^ 


(4-41) 


and  obtain  the  relation 


Ii  =  2kT  dxe 


-Ihi)  costf  x/kT  _  -1iw/2kT 

“  iCx  6 


MW  , 

coshy 


=  kTe“^/^*''^  K  (hu/2kT) 


(4-42) 


Thus  the  contribution  is  given  by 


<“>  •  [?  +  2  (^  +  ^)]  u e 


dipole 


(1 ■'ofe) 


(4-43) 
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We  note  further  that  since  llm^  ~  ttu/kT<^  o  that 


dipole 


(u)  vanishes  as  (tP,  or 


t  /■  A 

O'!,,  ,  (w) 

UH»o  dipole 


-  ^  iw ""  \\  “  ^  ^  ^  ^->]  ^ 


The  ratio  of  this  value  to  the  limit  for  small  hu/kT  of  the  static  contribution  Is 

(4-45) 


^^-Thr)  L-^(5^ - 


11m 

'‘static 
kT 


As  an  example,  consider  a  closed-shell  atom,  for  which  3^  -  =  0,  Then 

O  A  Q  1 

f'^lO  cm  (Flrsov  and  Chlblsov,  1960),  while  usually  0(0) ~  10  cm*.  Thus 
even  for  the  hl^^est  frequencies  we  might  consider,  the  dipole  contribution  to 
free-free  absorption  may  be  neglected  with  respect  to  the  static  potential  effect. 

Since  the  Bremsstrahlung  cross-section  for  the  molecule  had  the  same  form 
as  that  for  the  neutral  atom  In  our  approximation,  we  may  write  Immediately  for 
the  free-free  absorption  coefficient  In  the  field  of  molecules 


64  /Fn  n  aor(0)'h  (kT) 

^  ft  ft 


1/2 


molecule 


3u? 


m 


JJT 


(4-46) 


84 


THE  UNIVERSITY  OF  MICHIGAN 

4386-2-F 


4.10  Relation  to  the  Loreatz  Formulation 

The  interaction  between  radiation  and  an  ionized  gas  la  i^nomenologically 
described  by  the  coefficients  of  absorption  and  emission,  or  alternately  by  the 
electrical  c(mductivity.  The  alternate  use  of  these  coefficients  tends  to  obscure 
the  fact  that  all  of  them  are  manifestations  of  the  same  atomic  mechanisms.  The 
mechanism  pre-eminent  in  the  radio-frequency  range  is  free-free  transitions. 
(Paraphrase  of  Introductory  remarks  of  Oster,  1961) . 

This  discussion  refers  to  the  situation  of  operation  at  frequencies  above 
the  plasma  frequency,  where  at  least  the  classical  description,  upon  which  most 
engineers'  intuition  is  based,  attributes  attenuation  to  the  dissipative  effect  of 
collisions.  That  is,  the  Lorentz  attenuation  coefficient  (i.  e.  the  H  =  0  case  of  the 
Apideton-Hartree  equations)  vanishes  for  vanishing  collision  frequency.  R  is 
clear  that  our  quantum-mechanical  approach  is  describing  the  same  process  as 
the  collisional  dissipation  in  the  Lorentz  treatment.  Their  numerical  predictions 
should  then  agree  well  to  the  accuracy  of  the  approximations  inherent  in  both 
procedures,  and  their  dependence  on  u,  T,  n  and  should  be  closely  related. 

The  good  numerical  agreement  has  been  verified  in  a  number  of  calculations 
which  were  made  on  classified  problems  and  submitted  to  ARP  A  and  AFCRL,  the  spon¬ 
soring  agencies  for  this  contract.  In  order  to  preserve  the  unclassified  status  of 
thiw  report,  we  shall  not  reproduce  them  here.  It  is  of  interest,  however,  to  compare 
the  functional  form  of  the  attenuations  in  the  two  approaches. 
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After  some  manipulation^  the  absorption  coefficient  for  free-free  transitions 
in  the  static  field  of  neutral  atoms  (equation  4-38)  may  be  written 


ol.  (u)  = 


4  A 

c  ^  U)  '  * 


(4-47) 


where  we  have  defined 
4jrn  e* 

duoed  uf  - -  . 

p  m 


c  a\,  rmj 


0(0)  as  the  collision  frequency,  and  intro- 


Correspondingly,  the  power  attenuation  coefficient  from  the  Lorentz  ex¬ 
pression  is  approximated  to  order(— ^  by  the  expression 

"l*"*  '~t 

t  ''o\» 

when^-jj-y  <<lalso.  Now  we  expect  that  collisional  dissipation  or  free-free 


(4-48) 


absorption  should  be  the  predominant  process  only  a  bit  above  the  plasma  frequency, 
since  a  resonance  may  exist  there.  Thus,  it  is  reasonable  to  use  the  two  express-  ' 
ions  only  for,  rou^ily,  u J^.<  1/3.  Since  the  next  term  in  the  expansion  of  the 
Lorentz  equation  is  of  order  #  it  is  justified  to  use  the  approximate  form 

(4-48)  in  the  Joint  region  of  validity  of  both  approaches.  The  restriction  on  collision 
frequency  is  related  to  the  averaging  process  emjdoyed  in  the  derivation  of  the 
Lorentz  equations;  it  must  be  meaningful  to  speak  of  the  electron  as  between 
ooUlsicHis  during  a  cycle,  and  collisions  should  not  overlap  each  other.  An  un¬ 
certainty  principle  argument  which  demands  that  the  relative  uncertainty  in  photon 
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energy  be  « 1,  and  relates  this  to  uncertainties  in  electron  energy  corresponding 

to  measuring  this  within  a  At-  mean  free  time,  l/i/  ,  yields  the  same  restriction 

c 

on  for  the  quantum- mechanical  case.  We  thus  conclude  that  for  the  regimes 

in  which  both  expressions  may  safely  be  used,  they  have  the  same  functional 

4  /F 

dependence,  and  agree  to  within  a  factor  of  — ^ —  ^1.9.  This  factor  could  be  a 
result  of  imprcqper  averaging  in  the  Lorentz  derivation,  or  an  inaccuracy  of  the 
perturbation  theory  used  to  derive  the  quantum-mechanical  result.  In  any  case, 
both  formulas  are  clearly  talking  about  the  same  thing. 

For  the  case  of  electron-ion  interaction,  a  collision  frequency  to  be  inserted 
into  the  Lorentz  equation  is  given  by  (Sikarofsky,  1961) .  It  is,  in  our  notation. 


4 

3 


(4-49) 


Computations  using  this  eiqpression  agree  well  with  those  using  our  relaticm  (4-34). 
The  functional  dependence  is  the  same,  except  for  the  In  here  replacing  our  sum 
of  terms  inside  the  b  racket.  For  parameters  we  have  used,  these  correction 
factors  differ  little. 

Attenuation  processes  below  and  in  the  immediate  neighborhood  of  the 
plasma  frequency  have  not  been  studied;  the  following  statements  are  thus 
conjectural.  Intuition  based  on  the  Lorentz  treatment  suggests  that  attenuaticm 
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below  the  plasma  frequency  Is  purely  reactive  and  non-zero  for  a  collisionless 
plasma,  and  not  linearly  dependent  on  collision  frequency  when  this  quantity  is 
introduced  (Ratcll^e,  1959).  Free-free  absorption  by  uncorrelated  events  is 
probably  not  the  correct  picture  here,  and  we  suspect  that  the  Lorentz  equations 
could  be  in  error,  too.  An  investigation  based  on  the  microfield  in  a  plasma,  as 
well  as  the  classical  or  quantum-mechanical  dynamics  of  the  collisions,  seems  to 


be  the  appropriate  approach. 
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V 

CONCLUSIONS 

An  InveBtigation,  both  theoretical  and  experimental,  has  been  initiated  into 
the  absorption  of  electromagnetic  radiation  by  a  plasma  at  frequencies  exceeding 
the  plasma  frequency. 

The  theoretical  work  has  concentrated  on  free-free  absorption  rather  than 
on  any  resonant  effects  that  may  occur  in  the  immediate  vicinity  of  the  plasma 
frequency.  Ejqpressions  have  been  obtained  for  the  free -free  absorption  coefficient 
in  the  presence  of  both  ions  and  neutral  particles.  These  esqiressions  are  evalu¬ 
ated  e^qplicltly  for  the  case  of  a  Maxwellian  electron  distribution,  and  given  in 
terms  of  integrals  over  the  distribution  function  for  non-equilibrium  situations.  In 
the  latter  case,  since  the  distribution  function  depends  in  principle  on  the  spectral 
density  of  radiation  present  and  the  absorption  coefficient,  the  solution  is  purely 
formal.  However,  it  may  be  used  in  situations  for  which  the  distribution  functicm 
is  primarily  determined  by  such  effects  as  thermal  gradients  or  potential  gradients 
rather  than  by  interaction  with  incident  radiation. 

The  e^qperlmental  investigation  has  obtained  data  for  the  interaction  of 
radiation  with  a  pure  electron  stream.  A  simple  physical  model  based  on  Lorentz 
theory  and  introducing  a  collision  frequency  with  the  walls  succeeds  only  partially 
in  accounting  for  the  data,  hi  particular,  the  absorption  of  radiation  by  the 
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electron  stream  in  the  vicinity  of  the  plasma  frequency  is  not  e^lained  satisfac¬ 
torily  as  yet.  Attempts  to  introduce  gases  such  as  helium  have  caused  cathode 
poisoning  due  to  impurities  evolving  in  the  system  at  a  rate  more  rapid  than  ex- 
pectedi  For  this  reason  the  data  obtained  on  helium  from  the  first  cathode  was 
not  sufficiently  reliable  to  warrant  analysis.  Time  did  not  permit  helium  runs. 

I 

hi  order  to  interpret  the  eiqierimental  results  it  is  necessary  to  know  the 
electron  distribution  within  the  tube.  A  considerable  effort  has  been  devoted  to 
this  problem^  with  the  results  to  date  presented  here. 

The  theoretical  Investigation  of  free-free  absorption  led  to  the  conclusion 
that  for  frequencies  reasonably  above  the  plasma  frequency,  the  classical  Lorentz 
approach  predicts  the  attenuation  with  satisfactory  accuracy,  or  more  precisely, 
agrees  well  with  the  free-free  absorption  coefficient.  It  was  hoped  that  the  eiqieri- 
mental  results  would  provide  a  check  on  the  theoretical  work,  to  permit  a  decision 
on  whether  the  theoretical  results  are  both  correct,  or  merely  consistent  but  in 
error  because  of  the  neglect  of  other  attenuation  mechanisms.  This  hope  has  not 
yet  been  completely  realized,  and  we  feel  that  the  experiments  required  should  be 
performed. 

We  also  have  serious  doubts  concerning  the  validity  of  the  Lorentz  equa¬ 
tions  below  the  plasma  frequency,  and  feel  strongly  that  e}q)erimental  and 
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theoretical  work  in  this  area  should  be  initiated,  as  should  a  study  of  the 
immediate  nei^borhood  of  the  plasnoa  frequency,  where  the  interplay  of  several 
phenomena  creates  a  complex  situation. 
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APPENDIX  A 

A  VARIATIONAL  DETERMINATION  OF  THE  PROPAGATION 
CONSTANT  IN  AN  INHOMOOENOUS  COAXIAL  LINE 

Thif  i4>p«ndix  is  concerned  with  the  propagation  characteristics  of  a  coaxial 
wave  guide  which  is  filled  with  a  dielectric  substance  v/boBe  permittivity  varies 
only  radially.  At  present  we  shall  discuss  an  averaging  procedure  which  will 
enable  us  to  take  into  account  the  variable  nature  of  the  dielectric  substance.  Ulti¬ 
mately,  it  is  hoped  that  a  perturbation  procedure  may  be  derived  ^ich  will  enable 
us  to  correct  this  averaging  method  in  order  to  give  a  more  detailed  account  of  the 
presence  of  the  dielectric. 

We  consider  first  a  coaxial  wave  guide  which  is  filled  with  a  dielectric  of 
permittivity  e  -  e(r),  and  a  constant  permeability  With  the  time  dependence 
assumed  to  be  e]g>(Jut),  the  components  of  the  electric  and  magnetic  fields  may  be 
derived  from  the  vector  potential  w  by  the  following  relations 

H  =  ju  V  X  (ef) 

and  E  =  1/^  V  X  V  X  (ex) 

If  we  take  r  to  depend  on  a  single  component,  that  is  x  =  i^  itz> 

c  x^  +  =  0  (A-1) 
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From  the  e^resBions  E  and  H  we  may  calculate  the  e]q)licit  forms  of  the 

components  in  terms  of  ir„,  that  is, 

z 


-T,  -  1  9  9  .9^ 

r  ~  r  ar  ar  ^  "  dzdr  ^ 


1  a^ 

P  r  aza0  2 

„  ir9  9^.ia2  ] 

*  r  ar  ar  ^  r  a^2  zj 


^  ^  'z  ' 

8z2  *  °  * 


ju  a 

H  = - -r  ^  *■-  * 

r  r  dP  ® 


and 

H^  =  0 

where  the  coordinates  r,  0,  z  are  the  customary  cylindrical  coordinates.  The 
above  fields  are  divergence  free,  that  is 

V  •  (Mq  H)  =  0 

and 

V  '  (€E)  =  0 
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Furthermore,  these  components  describe  an  £  type  of  field,  there  being  no  need 
of  an  H  type  of  field  since  we  shall  require  excitation  in  terms  of  the  principal 
mode. 


Since  we  shall  ultimately  be  concerned  with  an  excitation  idiich  is  axially 
symmetric,  that  is  Independent  of  the  angle  0,  the  field  components  \^ich  we  just 
described  assume  the  form 


eEr  = 


c  E. 


dzdr 

1  a 


r  -  C  If  “ 

r  dr  dr  z 


dz^  z  Z 


Hwi  =  -  ju  -  ex 

P  ar  2 


and 


Consider  now  a  coaxial  wave  guide  of  inner  radius  a  and  outer  radius  b 
with  perfectly  conducting  walls.  On  these  walls  we  have  the  boundaiy  conditions 
£2  z)  =  (b,  z)  =  0.  The  presence  of  the  variable  permittivity  e  =  e(r) 

complicates  the  formation  of  solutions  of  equation  (A-1),  from  two  points  of  view. 
That  is,  the  e  is  contained  both  in  the  differential  equation  as  well  as  in  the  bound¬ 
ary  condition,  hi  order  to  circumvent  this  difficulty,  we  note  that  we  can  write 
separated  solutions  of  (A-1)  in  the  form 
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=  63?)  (jOfz)  f(r). 

From  this  we  find  that  f(r)  satisfies  the  ordinairy  linear  differential  equation 

^  —  r  +  (HqJ  -a^)f  =  0  (A -2) 

r  dr  dr 

and  (A -2)  is  to  be  solved  subject  to  the  boundary  conditions  f(a)  =  f(b)  =  0. 

Now  we  are  in  a  position  to  describe  a  variation  principle  for  determining 
the  smallest  characteristic  value  a.  That  is,  we  have  an  effective  method  of 
finding  the  smallest  characteristic  value  or  of  equation  (A -2)  solved  subject  to  the 
boundary  conditions  which  are  imposed  by  the  problem,  hideed,  upon  mult4>lying 
equation  (A-2)  throu^  by  rf(r),  integrating  by  parts  and  noting  the  boundary  con¬ 
ditions,  we  get 

J  c^(r)  r  f^(r)  -  rr^(r)  dr 

(  rf^dr 
•'a 

and  for  the  class  of  admissable  functions  f,  a  is  above  a  minimum  a^.  will  be 
a  minimum  for  the  correct  function  f. 
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-  A 

4  e  dx 

X 


dx 


or 


3 

2 


1  .X* 


-46  dx  +  - 


1_ 

2 


X* 

e 


L 


(B-2) 


Substituting  (B-2)  into  (B-1)  one  obtains 


1  1 
2  X 


4  X* 


(B-3) 


In  the  next  step  one  obtains 


; 


X*  J 

e  dx 


1  x^ 
-  e 

X 


1_ 

x" 


+4  r.» 


*1 V 

8  x^ 


15  f  ^ 

8  J  X* 


e^dx 


(B-4) 


This  procedure  can  be  continued  as  many  times  as  necessary  whereby  each  time 
the  integrand  of  the  remainder  integral  on  the  right  hand  side  is  decreased  by 
the  factor  — ^  .  For  x  »  1  the  remainder  integral  becomes  negligible  con^ared 
with  the  integral  on  the  left  hand  side  after  a  few  steps. 
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